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1.  Introduction 

When  Ice  crystals  fall  through  the  0°C  isotherm  in  the  atmosphere,  they 
begin  to  melt.  In  this  Ice  melting  zone,  it  is  known  that  radar  echoes  show 
the  so-called  "bright  band"  (Ryde,  1946).  The  bright  band  is  believed  to  be 
due  to  microwave  back  scattering  from  melting  ice  crystals  that  are  undergoing 
simultaneous  microphysical  changes  and  interactions. 

The  bright  bands  are  generally  observed  from  about  100  to  400  m  below  the 
0°C  isotherm  (Austin  and  Bemis,  1950;  Hooper  and  Kippax,  1950),  with  the  ratio 
of  the  radar  reflectivity  in  the  center  of  the  bright  band  to  that  at  the  0°C 
isotherm  being  normally  between  15  and  30.  The  ratio  ranges  as  high  as  100  and 
as  low  as  2  or  3.  The  ratio  of  the  radar  reflectivity  in  the  brightest  position 
of  the  band  to  that  in  the  rain  just  below  varies  normally  in  the  range  from  4 
to  9. 

The  main  microphysical  factors  that  are  thouqht  to  affect  the  charac¬ 
teristics  of  the  bright  bands  are  (1)  high  radar  reflectivity  of  water  and  its 
coating  effect  on  ice,  (ii)  fall  velocity  variation,  normally  the  increase,  (iii) 
variation  of  aggregation  process  Including  possible  disintegration  and  (iv) 
particle  shape  change  particularly  under  partially  melting  condition  (Battan, 
1973). 

The  laboratory  studies  of  ice  melting  In  the  past  were  mostly  with  spherical 
Ice  particles  (Mason,  1956;  Macklin,  1963;  Drake  and  Mason,  1966;  Bailey  and 
Macklln,  1%'j)  although  different  levels  of  Ice  density  were  sometimes  chosen 
(Mason,  1956)  and  snow  pellets  were  recently  used  (Matsuo  and  Sasyo;  19B2). 

Matsuo  and  Sasyo  (1981a)  also  reported  melting  of  snowflakes  that  were  supported 
on  nylon  nets  and  it  Is  the  only  work  of  the  kind  known  to  date.  The  basic  micro¬ 
physical  Informations  of  Ice  crystal  and  snowflake  melting  are  essential  for 
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clarification  of  the  bright  bawl  formation  mechanism  but  they  aise  presently  far 
from  sufficient  for  the  purpose. 

In  our  laboratory,  an  effort  has  recently  been  made  to  determine  some  of 
the  needed  microphysical  data  under  simulated  atmospheric  conditions  by  devel¬ 
oping  and  utilizing  a  new  vertical  wind  tunnel.  Applying  the  obtained  data  in 
the  tunnel  both  to  the  existing  and  newly  developed  theories  of  ice  crystal  and 
snowflake  melting  under  possible  atmospheric  conditions,  the  bright  band  for¬ 
mation  mechanism  was  reexamined.  Jrfe  shall  report  results  of  the  study  below. 


2.  Experimental  Method 

In  planning  the  measurement  of  ice  crystal  and  snowflake  melting,  the  degree 
of  simulation  to  the  natural  melting  layer  was  considered.  Natural  snow  crystals 
and  snowflakes  melt  while  they  are  falling  and  due  to  melting  the  fall  velocity 
changes.  The  velocity  change  in  turn  affects  the  melting  rate.  As  discussed 
in  the  introduction,  the  fall  behavior  of  melting  hydrometeors  is  already  known 
to  be  one  of  the  main  factors  that  influence  the  bright  band  characteristics. 

For  this  reason  the  vertical  wind  tunnel  method  was  selected  for  measurements. 

2.1  Vertical  wind  tunnel  and  the  supporting  devices 

Vertical  wind  tunnels  were  originally  developed  to  test  aircraft  spin 
characteristics  using  models.  They  are  featured  with  (i)  a  saucer-shaped  velocity 
profile  with  low  turbulence  In  the  test  section  which  is  about  5  -  10%  lower 
at  the  centerline  than  at  the  side  to  keep  the  model  near  the  center  of  the  tunnel, 
and  (11)  a  fast  response  drive  system  to  rapidly  adjust  the  air  velocity.  Re¬ 
searchers  have  used  various  modifications  of  the  vertical  wind  tunnel  in  cloud- 
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physical  studies,  drop  suspension  In  particular  (Klnzer  and  Gunn,  1951;  Telford 
et  al.,  1955;  Klnzer  and  Cobb,  1956  and  1958;  Cotton  and  Gokhale,  1967;  Irlbarne 
and  Klemes,  1970).  Hoffer  and  Mallen  (1968)  developed  a  tunnel  which  was  capable 
of  suspending  droplets  smaller  than  500  in  in  diameter  by  Incorporating  a  veloc¬ 
ity  profile  creator  and  a  diffuser.  The  most  widely  used  vertical  wind  tunnel 
to  date  Is  that  In  the  University  of  California  at  Los  Angeles  (UCLA)  (Pruppacher 
and  Nelburger,  1968;  Beard  and  Pruppacher,  1969;  Pruppacher  and  Beard,  1975; 
Pruppacher  and  Rasmussen,  1979).  Ice  particles  were  also  studied  in  the  wind 
tunnel  with  (List,  1959)  or  without  (Pitter  and  Pruppacher,  1973;  Pflaum  and 
Pruppacher,  1979)  support.  Freely  suspended,  however,  were  only  dense  ice  par¬ 
ticles  like  frozen  drops  and  graupel. 

Suspension  of  snow  crystals  In  the  tunnel  presents  additional  problems  due 
to  their  non-spherical  shapes.  Indeed,  Pruppacher  and  Neiburger  (1968)  felt 
that  aerodynamic  characteristics  of  a  melting  snow  crystal  or  a  snowflake  may 
make  it  Impracticable  to  freely  suspend  them.  Fukuta  et  al.  (1979)  were  among 
the  first  to  freely  suspend  ice  crystals,  though  small.  The  group  later  devel¬ 
oped  an  upward  converging  tunnel  to  suspend  larger  crystals  for  longer  periods 
of  time  (Kowa,  1981;  Fukuta  et  al.,  1982). 

The  design  of  the  tunnel  was  decided  considering  (1)  the  available  and 
proven  technique  of  the  vertical  wind  tunnel,  (11)  the  snow  crystal  fall  velocity 
and  the  range,  (111)  operator  reaction  for  air  speed  adjustment,  (1v)  the  mech¬ 
anism  to  Introduce  the  crystals  Into  the  tunnel  and  Its  feasibility  in  the 
selected  design,  and  (v)  data  collection  under  outdoor  operation.  The  vertical 
wind  tunnel  section  consists  of  the  vertical  tunnel,  the  settling  chamber,  the 
honeycomb,  screens,  and  the  air  exit  box  with  snow  Introducing  device.  (See 


Major  sections  of  the  experimental  apparatus 


Vertical  Wind  Tunnel 


A  vertical  wind  tunnel  employing  an  upward  diverging  design  was  built 
and  tested  independently  for  optimization  using  suspension  of  model  snow  crystals 
as  the  tunnel  performance  indicator.  The  model  snow  crystals  were  formed  by 
placing  a  sheet  of  heavy  duty  aluminum  foil  over  a  piece  of  chamois  and  punching 
out  disks,  the  result  is  a  5  mm  diameter  aluminum  disk,  curved  upward  slightly 
at  the  edge  and  weighing  about  2  mg. 

The  diverging  design  would  hold  a  velocity  profile,  increasing  towards  the 
bottom,  and  would  allow  the  snowflake  to  automatically  find  a  balancing  posi¬ 
tion  in  the  tunnel.  Then,  as  it  melted,  and  its  fall  velocity  increased  due 
to  its  smaller  profile,  it  would  slowly  settle  toward  the  bottom  of  the  chamber. 
Depending  on  its  size,  an  adjustment  in  the  chamber  airflow  may  not  be  neces¬ 
sary.  Required  adjustment,  if  any,  would  be  small.  In  this  way,  the  design 
of  the  tunnel  would  simplify  the  process  of  adjusting  the  airflow  to  hold 
the  snowflake  motionless.  After  testing  different  shapes  including  square 
cross  section,  it  was  found  that  a  vertical  wind  tunnel  with  a  round  cross 
section  and  dimensions  60  cm  high,  9.5  cm  in  diameter  at  the  bottom  and  13.6  cm 
at  the  top,  giving  a  taper  of  less  than  3  degrees  in  semi-vertical  angle, 
was  satisfactory  for  our  purpose.  The  chamber  was  constructed  by  rolling  a 
24  mil  cellulose  acetate  sheet  and  embedding  at  both  ends  in  6.5  mm  thick 
acrylic  plates.  In  between  the  top  and  bottom  plates  were  two  3.2  mm  thick 
acrylic  rings  which  provided  support  approximately  10  cm  from  each  end  (see 
Fig.  2). 

Settling  Chamber 

A  settling  chamber  acts  to  constrict  the  flow  through  the  tunnel,  while 
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Fig.  2.  Vertical  wind  tunnel  and  supporting  equipments. 


damping  out  some  large  scale  turbulence.  The  settling  chamber  actually  used 
measures  22  x  30.5  cm2  and  is  8  cm  deep,  and  Is  made  of  21  gauge  galvanized 
sheet  metal.  However,  it  was  found  that  disturbance  in  the  settling  chamber 
would  not  affect  the  tunnel  performance  due  to  excellent  flow  control  by  the 
honeycomb  we  adopted. 


Honeycomb 

After  testing  different  honeycombs  including  Hexcel  (Hexcel  Corp.,  Dublin, 
Ca.),  properly  packed  plastic  straws  of  suitable  diameter  were  found  to  be 
convenient,  inexpensive,  available  and  effective.  It  is  said  that  in  order 
to  obtain  a  laminar  flow,  the  length  to  diameter  ratio  of  10  is  sufficient. 

Two  new  methods  were  developed  with  respect  to  the  honeycomb  application  in 
the  tunnel.  The  first  is  the  scheme  to  produce  the  desired  air  velocity  pro¬ 
file.  The  scheme  of  using  straws  of  different  diameters  reported  by  Hoffer 
and  Mallen  (1967)  was  found  to  give  too  rapid  velocity  variation  in  the  lateral 
direction.  Against  this  problem,  we  developed  a  scheme  of  varying  the  length 
of  straws.  Since  straws  give  resistance  to  the  air  flow,  the  longer  the  starws, 
the  slower  the  air  velocity  after  passing  through  them.  To  obtain  the  air 
velocity  profile  gently  concave  upward,  the  length  of  the  straws  has  to  be 
larger  at  the  center  (see  Fig.  3).  The  second  is  the  method  to  cut  straws 
easily  and  after  cutting  to  keep  them  together.  An  electrically  heated  wire 
mounted  on  a  piece  of  plywood  was  pushed  through  the  straws,  melting  them 
flush  with  the  top  of  the  honeycomb  section  and  sealing  them  together.  This 
allowed  for  easy  removal  of  the  honeycomb.  The  honeycomb  used  for  measurement 
is  9.5  cm  in  diameter,  made  of  3.2  mm  diameter  plastic  straws,  arranged  in  12 
levels,  with  the  straws  verylng  In  length  from  6  to  9.4  cm.  The  straws  were 
arranged  in  place  in  a  9.5  cm  inside  diameter  (ID)  clear  acrylic  tube,  10  cm 


Honeycombs  of  different  design  (a)  and  a  honeycomb  positioned 
in  the  acrylic  cylinder  (b). 
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long,  which  was  held  at  each  end  by  a  30  cm  square  piece  of  6.4  mm  acrylic 
sheeting  (see  Fig.  3). 

Screens 

While  the  honeycomb  was  useful  In  approximating  laminar  flow,  further 
refinement  of  the  air  stream  was  required  to  achieve  the  desired  velocity  pro¬ 
file.  It  was  done  with  screens  applying  the  trial  and  error  method. 

According  to  Pankhurst  and  Holder  (1952),  screens  act  to  remove  large- 
scale  eddies  at  the  expense  of  the  Introduction  of  a  larger  nwnber  of  smaller 
eddies  which  decay  rapidly.  The  screen  thus  decreases  the  turbulence  at  a 
sufficient  distance  downstream,  although  it  may  considerably  increase  it 
immediately  after  the  screen.  Pope  and  Harper  (1966)  state  that  screens  are 
far  more  useful  than  honeycombs  for  adjusting  airflow.  Low  turbulence  wind 
tunnels  may  use  six  or  more  whereas  a  regular  tunnel  will  only  incorporate  one 
or  two.  According  to  Pope  (1954),  the  spacing  between  screens  doesn't  seem 
to  make  any  difference  as  long  as  they  don't  touch  each  other.  Our  screens 
are  made  using  18  x  16  aluminum  mesh,  mounted  on  rings  constructed  of  20  gauge 
steel  wire. 

Three  layers  of  screens  were  used,  directly  above  the  honeycomb  in  our 
tunnel.  The  desired  velocity  profile  was  made  by  using  screens  of  higher 
reslstence  in  the  center.  These  screens  were  prepared  by  marking  a  circular 
center  section  on  the  mesh  and  then  removing  the  excess  wire  strands.  The  re¬ 
maining  strands  were  arranged  In  the  radial  direction,  evenly  spaced  and  attached 
to  the  ring  with  epoxy  glue.  The  rings  and  mesh  were  spray  painted  black  and 
the  rings  further  coated  with  polyurethane  paint.  A  velocity  profile  whose 
center  dip  was  too  shallow  allowed  the  artificial  snow  crystal  to  work  its  way 
up  the  slide  until  gravity  pulled  It  down.  It  would  cross  the  bottom  of  this 


velocity  MdishH  and  climb  higher  on  the  opposite  side.  By  continuing  this  pro¬ 
cess,  the  artificial  snowflake  would  soon  climb  over  the  side  of  the  velocity 
"dish"  and  Impact  with  the  suspension  chamber  wall.  If  the  velocity  profile 
had  been  too  deep,  the  artificial  snow  crystal  drifted  out  of  the  bottom  of 
the  "dish"  and  tumbled  Instead  of  simply  sliding  back.  Once  the  crystal  began 
tumbling.  It  was  Impossible  to  hold  It  In  the  middle.  Fig.  4  shows  the  screens 
used  and  their  order  of  placement  on  the  suspension  chamber,  beginning  with  the 
upper  left.  The  height  of  the  screens  above  the  bottom  of  the  tunnel  is  as 
follows: 


first 

0 

second 

3.6 

third 

8.0 

fourth 

10.6 

fifth 

13.0 

sixth 

14.1 

seventh 

15.9 

The  best  results  were  achieved  by  placing,  bottom  to  top,  the  middle,  small 
then  large  screen,  and  not  the  small,  middle  and  large  as  one  might  expect. 

Snow  Introducing  Device 

To  receive  snow  crystals  during  windy  conditions,  a  chimney  was  erected 
based  on  experiences  by  Zlkmunda  and  Vali  (1972).  It  was  a  cylinder  24  cm 
in  diameter,  91.4  cm  high  made  of  heavy  wire  mesh  and  covered  with  1.5  mil 
plastic  sheet.  This  cylinder  was  anchored  to  a  wooden  base,  and  positioned 
outside  the  shelter  and  above  the  vertical  wind  tunnel.  The  chimney  was  deep 
enough  to  enable  a  snowflake,  to  regain  Its  vertical  fall  in  the  tunnel  even 
If  It  had  previously  been  travelling  almost  horizontally,  thus  allowing 
measurements  on  windy  days.  On  the  top  of  this  tunnel,  opening  controls  (3 
removable  covers  with  center  holes  3  to  15  cm  in  diameter)  can  be  attached. 
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Depending  on  the  intensity  of  the  snowfall,  these  allow  finer  control  on  the 
number  of  snowflakes  entering  the  tunnel. 

The  snowflake,  after  falling  through  the  chimney,  reaches  the  top  of  the 
tunnel,  where  the  air  exit  chamber  is  located.  This  is  an  area  of  great  tur¬ 
bulence,  which  could  cause  the  snowflake  to  be  thrown  against  the  chamber  wall 
and  lost.  To  help  the  snowflake  fall  past  this  area,  an  insertion  tube  was 
designed  and  constructed.  This  consists  of  a  6.4  cm  ID  clear  acrylic  tube 
which  extends  23  cm  down  into  the  tunnel  space.  This  allows  the  snowflake  to 
pass  through  the  high  turbulence  area  undisturbed.  There  is  some  turbulence 
at  the  bottom  of  this  tube,  but  not  enough  to  adversely  affect  the  snowflake 
before  it  falls  into  the  velocity  profile  and  attains  stability.  The  entrance 
to  this  insertion  tube  is  a  2.5  cm  diameter  opening  which  acts  to  center  the 
snowflake  in  the  tube.  The  top  opening  has  a  spring-loaded  cover  which  is 
opened  by  pulling  a  cable.  As  soon  as  a  snowflake  has  entered  the  tunnel  and 
achieved  stability,  the  cable  is  released  and  the  entrance  closes. 

Air  Exit  Box 

An  exit  was  designed  and  constructed  for  the  air  to  exit  smoothly  from  the 
tunnel.  It  is  30  cm  square,  11  cm  high  and  is  connected  with  the  valve  by  a 
4.1  cm  ID  clear  acrylic  tube.  Inside  this  box  is  a  15  cm  diameter  clear 
acrylic  cylinder,  which  is  sealed  to  the  top  and  bottom  of  the  box.  The  bottom 
of  the  box  is  cut  away  inside  this  cylinder  and  the  top  has  a  rubber  seal 
through  which  the  insertion  tube  fits.  This  15  cm  diameter  cylinder  has  500 
3.8  mm  diameter  holes,  spaced  6  mm  apart,  perforating  it.  To  insure  a  symmet¬ 
rical  airflow,  this  15  cm  diameter  cylinder  is  wrapped  with  a  1.3  cm  thick  foam 
rubber.  This  provides  a  resistance  to  the  airflow,  creating  an  area  of  lower 
pressure  in  the  box,  which  in  turn  causes  the  air  to  be  drawn  from  the  suspension 
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chamber  symmetrically,  rather  than  favoring  the  area  where  the  tube  to  the 
valve  is  located.  This  symmetrical  exit  of  the  air,  combined  with  the  center 
"dead  spot"  form  the  insertion  tube  tends  to  enhance  the  concave  upward 
velocity  profile  in  the  tunnel  (see  Fig.  5). 

2.2  Air  suction  system 

Air  Suction 

The  air  suction  for  our  apparatus  is  provided  by  16  gallon  (60  liter) 

Wet-Dry  Shop  Vac  (Sears,  Roebuck  Co.).  It  provides  more  than  adequate  flow 
rates  for  this  study  and,  is  rugged  enough  to  withstand  cold  temperatures  for 
long  periods  of  time.  In  addition,  it  met  the  portability  requirements  of  our 
study  and  attachments  and  replacement  parts  are  readily  available.  A  variac 
was  used  to  reduce  the  suction  rate  so  that  the  airflow  in  the  tunnel  fell 
within  the  range  suitable  for  suspending  snowflakes.  The  power  for  our  entire 
apparatus  was  provided  by  a  portable  alternator  (Sears,  Roebuck  and  Company, 

2600  W).  This  unit  provided  ample  power  to  meet  our  needs.  It  was  easily 
handled  by  two  people  and  operated  well  in  the  adverse  weather  conditions  in 
which  we  carried  out  our  measurements. 

Valve 

Direct  control  of  the  tunnel  air  velocity  is  provided  by  a  valve  designed 
and  constructed  In  our  laboratory  (see  Fig.  6).  This  valve  Is  located  directly 
outside  the  exit  portion  of  the  tunnel,  thereby  minimizing  the  lag  time  between 
valve  adjustment  and  the  corresponding  air  velocity  change  in  the  tunnel.  The 
valve  Is  constructed  of  a  10.8  cm  section  of  clear  acrylic  tubing  with  ID  of 
14.6  cm  and  outside  diameter  (OD)  15.2  cm.  The  Intake,  i.e.  exit  for  the  tunnel. 
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Fig.  5.  Air  exit  box  at  the  top  of  the  vertical  wind  tunnel. 


Fig.  6.  Valve  with  anemometer  probe. 


is  made  of  a  10  cm  section  of  clear  acrylic  tubing  with  5.7  cm  00.  The  bypass 
tubing,  coming  from  the  settling  chamber,  is  a  section  of  6.4  cm  diameter  vacuum 
hose  and  connects  to  the  valve  next  to  the  intake  previously  mentioned.  The 
actual  openings  into  the  valve  for  both  of  these  tubes  are  holes  3.3  cm  in  diam¬ 
eter.  These  holes  are  arranged  so  that  as  the  cover  plate,  made  of  polyethy¬ 
lene  to  slide  easily,  inside  the  valve  is  rotated,  the  amount  of  air  allowed 
to  enter  the  valve  remains  constant,  while  one  is  opening  the  other  closing. 

This  ensures  that  Shop  Vac  Is  drawing  a  constant  air  flow  rate  at  all  times, 
and  the  flow  rate  of  air  drawn  through  the  air  conditioning  system  also  remains 
constant  which  avoids  variations  in  the  relative  humidity  and  temperature  of 
the  air.  A  steel  shaft,  attached  to  the  cover  plate,  extends  through  the  top 
of  the  valve  and  outside  the  enclosing  box  and  is  attached  to  the  handle, 

56  cm  long.  The  valve  was  designed  so  that  an  upward  movement  of  the  handle 
results  in  increased  airflow  and  a  downward  movement  in  reduced  air  flow, 
making  the  tunnel  operation  easier. 

2.3  Air  conditioning  system 

The  original  design  of  our  apparatus  called  for  the  air  to  be  recirculated, 
giving  precise  control  over  the  relative  humidity  and  temperature.  Upon  actu¬ 
ally  testing  this  system  though,  two  serious  problems  appeared:  vibration  and 
heating,  both  caused  by  the  suction  system.  So,  it  was  decided  to  adapt  an 
open  system  passing  the  air  through  the  entire  system  one  time  only.  For¬ 
tunately,  the  temperature  of  the  environmental  air  remained  constant  over  the 
time  of  observation,  and  under  the  condition  a  constant  heating  and  airflow 
rate  through  a  humidifier  provided  satisfactory  control  of  air  conditions  in 
the  tunnel . 
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Heater 

The  heater  measures  44  x  55  cm2  in  a  wooden  frame.  Thirty  gauge  alumel 
wire  with  a  resistence  of  11  ohms  per  meter  is  strung  in  the  frame.  The  actual 
area  spanned  by  the  wire  is  40  x  46  cm2.  To  avoid  heat  loss  from  the  heater 
by  convection,  a  cover  with  holes  whose  cross-sectional  areas  combined  to  equal 
that  at  the  bottom  of  the  suspension  chamber  was  placed  over  the  heater  (see 
Fip.  7).  Over  this  cover  Is  another  box  which  draws  outside  air  through  a  6.4 
cm  diameter  flexible  tube.  The  metal  walls  of  the  conditioning  system  were 
insulated  with  1.9  cm  thick  styrofoam  sheeting. 

Humidifier 

The  relative  humidity  of  the  air  entering  the  suspension  chamber  was  con¬ 
trolled  by  an  evaporative  cooler.  Once  the  air  sample  has  passed  through  the 
heater,  it  can  be  directed  into  the  humidifier,  bypass  this  and  go  directly  to 
the  settling  chamber  or  any  combination  of  these  two.  This  allows  the  study 
of  snowflake  melting  in.  different  humidity  conditions. 

The  humidifying  section  comprises  the  bulk  of  the  conditioning  system 
(see  Fig.  7),  which  measures  44  x  55  cm2  and  is  79  cm  high.  It  has  three  layers, 
each  one  composed  of  two  layers  of  evaporative  cooler  batting.  It  holds  3.3 
gallons  (12.5  liters)  of  water.  After  running  for  several  minutes  (allowing 
time  for  the  batting  to  become  saturated),  this  leaves  8.6  cm  of  water  standing 
in  the  bottom.  The  batting  Is  held  in  place  by  0.5  In.  (1.3  cm)  wire  mesh, 
supported  by  wooden  corner  posts  and  plexiglass  crossmembers. 

The  Irrigation  system  is  made  up  of  a  standard  115  V  evaporative  cooler 
pump,  20  cm  in  height.  The  water  is  distributed  through  an  irrigation  system 
comprised  of  2.5  m  of  6.4  mm  diameter  rubber  hose,  with  holes  0.8  mm  in 


conditioning  system. 


diameter,  spaced  every  3.8  cm.  This  system  spreads  the  water  evenly  over  the 
top  layer  of  batting  where  gravity  draws  it  back  to  the  pump,  wetting  the 
batting  in  between.  This  system  is  very  effective  in  helping  achieve  maximum 
practical  relative  humidity. 

Above  the  irrigation  system  is  a  metal  cover  slide  to  control  humidifying 
or  bypassing  the  incoming  air.  The  relative  humidity  of  outdoor  air  is  re¬ 
presentative  of  that  beneath  a  storm  cloud  just  before  or  at  the  beginning  of 
snowfall.  As  the  snowfall  continues,  melting  snow  evaporates  in  the  air, 
raising  the  relative  humidity.  Under  the  condition,  simultaneous  control  of 
the  metal  cover  and  the  bypass  valve  will  permit  the  relative  humidity  ad¬ 
justment  to  the  desired  level.  To  insure  thorough  mixing  of  this  "dry"  air 
with  air  passing  through  the  humidifier,  a  device  was  built  (see  Fig.  7). 

The  "dry"  air  exits  the  bypass  line  through  ten  2.5  cm  diameter  polyvinyl 
chloride  (PVC)  pipes,  each  38  cm  long  and  containing  thirty-six  3.2  mm  diam¬ 
eter  holes.  The  humidified  air  passes  around  these  pipes  prior  to  being  drawn 
into  the  settling  chamber  and  mixes  with  the  dry  air. 

2,4  Photographic  section 

The  photographic  system  used  in  the  present  study  is  shown  in  Figs.  8 
and  9.  Through  the  use  of  a  common  shutter  release,  both  the  cameras  are  acti 
vated  and  a  timing  mark  is  placed  on  the  recorder  chart  simultaneously.  This 
camera  arrangement  gives  a  top  and  side  view  of  the  snowflake  at  desired  in¬ 
tervals  throughout  the  melting  process.  The  cameras  are  Olympus  OM-2  35  mm 
cameras,  equipped  with  Olympus  Winder  2s.  To  provide  further  enlargement  of 
the  snowflake,  the  top  camera  uses  a  14  mm  extension  tube  and  the  bottom  a 
7  mm  one.  Ilford  HP5,  ASA  400  film  is  used  in  the  cameras.  The  background 
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8.  Photographic  section 


lighting  is  provided  by  a  quartz-halogen  automobile  driving  light.  Intermit¬ 
tent  lighting  is  obtained  by  an  Olympus  T32  electronic  flash,  connected  elec¬ 
trically  to  the  top  camera  and  positioned  as  shown  in  Figs.  8  and  9.  To 
reduce  glare  from  the  tunnel  wall,  arrangements  were  made  whenever  necessary. 

2.5  Supporting  equipment 

Recorder 

To  record  the  outside  air  temperature,  the  time  the  photographs  were  taken, 
the  temperature  and  the  velocity  of  the  air  in  the  wind  tunnel  a  three-channel 
recorder  (Cole-Parmer,  model  8373-30)  was  used.  By  connecting  the  camera 
system  to  a  temperature  recording  channel  in  such  a  way  as  to  short  circuit 
it  when  a  photograph  was  taken,  all  the  above  parameters  were  recorded  on  it. 

Anemometer 

The  probe  of  an  anemometer  (Kurz  Instruments  Inc.,  model  430  M-4)  was 
mounted  in  the  tube  which  led  from  the  exit  of  the  suspension  chamber  (Fig.  6), 
to  utilize  the  sensitive  zone  of  the  meter  for  calibration.  With  the  probe 
location  we  chose,  calibrations  were  made  to  six  horizontal  reference  lines 
in  the  wind  tunnel,  2.5  cm  apart,  beginning  18.4  cm  above  the  bottom  of  the 
chamber.  These  lines,  in  conjunction  with  a  photograph,  allowed  an  exact 
determination  of  the  fall  velocity.  The  error  induced  by  the  anemometer- 
recorder  combination  was  1.2%. 

Hygrometer 

The  relative  humidity  of  sample  air  was  made  using  a  dew  point  hygro¬ 
meter  designed  and  built  in  our  laboratory.  The  mirrored  container  holding  ether 
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was  cooled  by  evaporation  as  the  rubber  bulb  was  squeezed,  forcing  air  bubbling 
up  through  it.  When  the  ether  cooled  the  container  to  the  dewpoint  temperature, 
condensation  would  form  on  the  mirrored  surface.  This  was  very  noticeable  with 
the  aid  of  the  uncooled  mirrored  fins  on  either  side.  The  dewpoint  and  sample 
air  temperatures  were  measured  by  two  calibrated  thermometers. 


Air  Temperature 


The  temperatures,  both  outside  and  inside  the  apparatus,  were  measured 
by  two  pre-calibrated  thermocouples,  constructed  from  type  T,  20  gauge  copper 
constantan  wire.  One  thermocouple  was  located  in  the  wind  tunnel.  The  thermo¬ 
couple  measuring  the  outside  air  temperature  extended  through  the  shelter  wall 
approximately  76  cm  and  was  protected  from  falling  snowflakes  by  a  cover. 

Shelter 

A  suitable  shelter  was  necessary  to  protect  both  the  equipment  and  the 
operators  from  the  environment,  as  the  observations  were  taken  during  snow¬ 
storms,  often  at  night,  with  temperatures  below  0°C,  and  lasted  for  several 
hours  at  a  time.  A  shelter  was  made  by  building  a  wooden  frame  on  the  back 
of  a  1968  Dodge  4-wheel  drive  Powerwagon.  The  framework  was  covered  with  two 
layers  of  plastic  sheeting,  and  a  layer  of  canvas,  held  in  place  by  ropes. 

The  equipment  was  mounted  inside  and  access  holes  cut  in  the  plastic.  During 
experiments,  the  ropes  were  removed  and  the  canvas  folded  back  to  expose  the 
access  holes  (holes  to  allow  the  snowflake  entrance  to  the  chamber  and  an  air 
intake  for  the  air  conditioning  system).  It  offered  the  further  advantage  of 
being  mobile,  which  allowed  travel  up  into  the  mountains  when  snow  was  not 
expected  in  the  valley.  The  shelter  during  an  observation  is  shown  in  Fig.  10. 


the  field 


2.6  Experimental  procedures 
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Preparation 

Due  to  the  limited  snowfall  during  the  winter,  mobility  was  an  important 
factor.  With  the  eventual  placement  of  the  experimental  apparatus  in  the  bed 
of  a  4-wheel  drive  truck,  travel  to  the  mountains,  where  snowfall  frequentl j 
occurred,  was  possible.  The  majority  of  the  equipment  was  stored  on  the 
truck.  When  it  was  determined  that  snowfall  would  occur  in  the  mountains, 
a  check  list  was  used  to  insure  everything  required  to  collect  data  was  loaded. 

The  problem  of  static  electricity  on  the  tunnel  surface  was  noticed  early 
in  our  equipment  development.  This  was  solved  by  wiping  the  surface  with 
damp  "Cling  Free"  antistatic  cloth,  drying  for  about  15  to  20  minutes,  and 
polishing  it  with  dry,  clean  cotton  rag.  It  was  found  that  this  treatment  of 
the  chamber  was  not  permanently  effective.  Therefore,  we  treated  the  equipment 
for  static  electricity  prior  to  every  storm  and  checked  for  its  presence  at 
the  beginning  of  each  data  collection  period. 

Equipment  Set  Up 

lJ|jon  arrival  at  the  data  gathering  site,  the  truck  was  parked  in  the  most 
level  place  available.  After  unloading  the  equipment,  a  decision  was  made  as 
to  whether  or  not  to  use  the  opening  controls  and/or  entrance  tunnel.  This  was 
based  on  the  snowfall  intensity  and  experience.  The  equipment  worked  best  when 
about  four  to  ten  snowflakes  entered  the  chamber  every  minute.  With  higher 
snowflake  concentrations,  It  became  difficult  to  concentrate  on  a  single  flake. 
Also,  the  others  could  Interfere  with  Its  suspension  and  eventually  land  on 
the  chamber  walls  or  screens  where  they  left  water  deposits.  These  water 
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deposits  could  disturb  the  air  velocity  profile  or  the  photography,  depending 
on  their  position.  Once  the  equipment  was  unloaded,  the  work  was  divided  be¬ 
tween  an  outside  and  Inside  worker. 

An  Important  initial  step  was  the  tunnel  stability  test  using  the  arti¬ 
ficial  snowflake.  Experience  showed  that  if  the  artificial  snowflake  could 
not  be  suspended,  the  chance  of  suspending  natural  snowflakes  was  very  small. 
This  test  was  very  useful  in  identifying  potential  problems.  After  the  sus¬ 
pension  test,  the  snowflake  was  removed  by  suction. 

The  Inside  worker,  like  the  outside,  also  had  many  tasks  to  complete. 

The  two  workers  had  to  communicate  as  they  set  up  and  confirm  the  proper  oper¬ 
ation  of  the  equipment.  The  camera  with  appropriate  extension  tubes  were 
placed  and  made  ready.  The  focal  points  of  the  cameras  were  In  the  center  of 
the  chamber,  between  the  third  and  fourth  reference  lines.  The  recorder  channel 
were  zeroed.  The  checklists  were  reviewed  just  prior  to  data  collection  to 
confirm  that  all  tasks  had  been  completed. 

Data  Collection 

As  with  setting  up  the  apparatus,  data  collection  was  also  a  two-man  job, 
with  one  person  (hereafter  called  controller)  operating  the  wind  tunnel  valve 
and  the  cable  release/marker  device  and  the  other  (the  supporter)  the  opening 
to  the  insertion  tube  and  taking  notes. 

Before  data  collection  began,  the  controller  measured  the  air  and  dewpoint 
temperatures  of  the  tunnel  while  the  supporter  was  recording  relevant  data. 

Data  collection  began  with  the  supporter  opening  the  entrance  to  the  insertion 
tube.  When  a  snowflake  entered  the  suspension  chamber,  he  Immediately  released 
the  control  cable,  whereupon  the  spring  loaded  opening  mechanism  closed.  Where¬ 
as  before  this  operation,  the  air  was  being  drawn  both  from  the  vertical  wind 


tunnel  and  from  the  outside,  now  It  came  solely  from  the  wind  tunnel.  This 
gave  an  Immediate  velocity  Increase  In  the  tunnel  which  was  enough  to  cause 
most  snowflakes  to  hesitate  momentarily  In  their  falling.  During  this  hesi¬ 
tation,  the  controller  Increased  the  velocity  and  gained  control  of  the  snow¬ 
flake  suspension.  When  the  snowflake  was  under  control,  approximately  two 
seconds  after  entering  the  tunnel,  the  controller  activated  the  camera  release/ 
marker  device  to  record  It.  He  also  increased  the  airflow  to  keep  the  snowflake 
vertically  stable  as  it  melted  and  Its  fall  velocity  increased.  This  required 
practice  to  keep  the  snowflake  In  the  optimum  position  for  photographing.  As 
the  melting  process  progressed,  both  workers  carefully  observed  the  behavior  of 
the  snowflake.  When  the  suspension  was  lost,  the  supporter  recorded  the  in¬ 
formation  concerning  that  individual  snowflake  on  the  recorder  strip.  The 
number  of  pictures  taken  was  also  recorded.  Then,  the  same  procedure  was  re¬ 
peated.  Cleaning  of  the  tunnel  was  carried  out  whenever  it  became  necessary 
following  the  predetermined  procedure. 

Throughout  our  data  collection  the  temperature  inside  the  chamber  ranged 
from  1.8  to  3.9°C,  without  operating  the  heater,  due  to  the  outside  air  tem¬ 
perature,  body  heat,  the  electric  heater,  and  other  equipment  inside  the  shelter. 
This  was  rather  convenient  as  these  were  some  of  the  temperatures  at  which  we 
desired  to  collect  data.  For  other  temperatures,  the  supporter  must  also  con¬ 
trol  the  heater  through  the  on-off  switch  located  on  the  front  of  the  instrument 
box.  By  reading  the  temperature  on  the  recorder  and  controlling  the  heater 
through  Its  switch,  temperature  could  be  maintained  at  any  level  within  ±0.2°C. 

Shutting  Down 

The  procedures  for  repacking  the  equipment  were  also  divided  between  the 
inside  and  the  outside  worker,  using  a  checklist  prepared. 
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Upon  returning  to  the  university,  the  supply  levels  were  checked.  If  data 
collection  was  expected  within  the  next  24  hours,  the  wind  tunnel  was  also 
serviced  for  static  electricity  as  described  earlier. 

Post  Experimental  Maintenance 

Thorough  cleaning  of  the  tunnel  is  crucial  for  the  successful  suspension 
of  ice  crystals.  Unfortunately,  cleaning  of  the  equipment  required  the  removal 
of  the  wind  tunnel,  the  honeycomb  device,  and  the  top  filter  section.  After  the 
device  was  removed,  water  sitting  in  the  tunnel  was  wiped  off  with  a  sponge  or 
rag.  Compressed  air  was  used  to  remove  any  water  that  had  settled  inside  the 
honeycomb  device.  The  final  procedure  involved  polishing  the  wind  tunnel  with 
the  antistatic. 

2.7  General  remarks  of  measurement 

The  outdoor  measurements  of  natural  snow  crystals  provided  us  with  some 
unusual  and  interesting  experiences.  It  may  be  worthwhile  to  mention  some  of 
them  here. 

For  thermocouple  operation  of  temperature  measurement,  an  ice  bath  of  0°C 
is  commonly  used.  Crushed  ice  will  normally  melt  to  yield  the  condition.  Under 
our  operating  conditions,  though,  the  crushed  ice  was  found  sometimes  not  to 
melt  due  to  low  air  temperatures.  So,  we  had  to  flush  the  crushed  ice  with 
water  and  keep  some  water  with  it. 

The  position  of  a  snowflake  in  the  wind  tunnel  was  important.  Unless 
centered  after  entering  the  upper  portion  of  the  tunnel,  they  hit  the  wall. 

In  the  upper  regions  of  the  tunnel,  the  air  velocity  profile  was  flattened  and 
the  snowflake  could  impact  with  the  wall.  When  centered,  the  velocity  profile 

moved  it  away  from  the  wall  as  it  fell  into  the  observing  area  and  was  suspended. 
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Snow  crystals  could  be  suspended  almost  anywhere  in  the  tunnel  depending  upon 
their  size.  Very  small  ones,  approximately  one  millimeter  or  less  in  diameter 
and  appearing  only  as  a  sparkle  of  light,  could  even  be  suspended  on  the  crest 
of  the  velocity  profile. 

It  was  found  that  our  apparatus  was  not  suitable  for  all  types  of  ice 
particles.  Graupel  proved  to  be  too  dense  to  suspend  and  snow  which  was  mostly 
melted  was  the  same.  Half  melted  snow  often  left  water  deposits  which  inter¬ 
fered  with  the  tunnel  airflow  and  our  visibility.  Best  results  were  obtained 
by  waiting  until  later  in  the  storm  when  the  surface  temperature  dropped  and 
the  snowflakes  were  not  melting.  Also,  moderate  snowfall  was  required  to  receive 
a  significant  number  of  snowflakes  into  the  tunnel. 

Large  snowflakes,  greater  than  ten  millimeters  in  diameter,  were  seldom 
observed  in  the  tunnel  due  to  several  reasons:  The  concentration  of  these  size 
flakes  during  a  storm  was  relatively  low.  Also,  due  to  small  entrance  to  the 
tunnel,  only  2.5  cm  in  diameter,  large  snowflakes  often  struck  the  side  of  the 
entrance  opening  causing  showers  of  small  crystals.  Once  in  the  tunnel,  the 
positioning  of  the  snowflake  became  a  factor. 


3.  Results  of  Measurement 

In  this  section,  we  describe  and  discuss  the  results  of  our  measurement, 
starting  from  our  visual  observation  during  the  tunnel  operation  in  the  field. 

3.1  Observation  of  melting  behavior 

The  data  collection  for  this  study  was  performed  during  the  three  storms 
In  January,  1982.  During  this  time,  nearly  400  photographs  of  Ice  crystals 
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and  snowflakes  were  taken  and  over  100  sequences  of  their  melting  were  analyzed. 
In  addition,  several  experiments  were  conducted  and  ice  crystal  melting  was 
observed  although  equipment  malfunction  prevented  photographic  data  collection. 

In  all,  about  2,700  ice  crystals  and  snowflakes  were  suspended  and  observed  to 
melt.  The  data  for  single  crystals  were  analyzed  and  reported  by  Savage  (1982) 
and  those  for  snowflakes  by  Donovan  (1982).  From  these  observations,  several 
interesting  patterns  surfaced. 

Rotation 

Usually  only  in  small  crystals,  from  1-3  mm  in  diameter,  we  observed  that 
the  crystals  rotated  around  their  vertical  axes.  He  called  this  "helicoptering". 
As  time  progressed,  the  crystal  began  orbiting  around  the  center  of  the  chamber, 
rather  than  helicoptering.  With  each  pass,  the  circle  became  wider  and  wider, 
until,  just  fractions  of  a  second  before  melting,  the  crystal  impacted  with  the 
chamber  wall.  This  spiral  was  observed  only  at  the  end  of  the  melting  process, 
and  only  with  this  size  crystal.  The  size  of  smaller  crystals  and  their  appear¬ 
ance  as  a  sparkle  of  light  prevented  us  from  observing  the  helicoptering  but 
larger  ones  did  not  show  the  rotation.  The  helicoptering  movement  did  not 
affect  our  ability  to  suspend  a  crystal.  Helicoptering  was  not  observed  in 
snowflakes. 

Vertical  Stability 

The  vertical  stability  of  a  suspended  snowflake  or  an  ice  crystal  was  de¬ 
pendent  upon  its  size  and  shape.  The  small  crystals,  less  than  one  millimeter 
in  diameter,  were  the  easist  to  suspend.  These  weren't  dependent  upon  the 
velocity  profile  to  hold  them.  They  could  be  suspended  anywhere  in  the  chamber 
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for  one  or  two  seconds,  during  which  time  they  melted.  Snowflakes  and  crystals 
became  increasingly  unstable  with  size,  due  to  the  way  they  melted.  Also  in 
our  wind  tunnel,  we  were  able  to  suspend  only  the  aerodynamically  stable  crys¬ 
tals.  Although  most  of  the  ice  crystals  which  entered  in  the  wind  tunnel  appeared 
to  be  aerodynamically  stable,  we  were  able  to  suspend  only  roughly  10*  of  the 
crystals  because  the  rest  of  them  fell  into  the  'instable  zones.  It  still  took 
much  skill  and  luck  to  suspend  even  those  crystals  which  fell  within  the  stable 
zones  of  the  tunnel . 

Melting  and  Disruption 

The  surprising  observation  during  this  study  was  the  way  snowflakes  and 
crystals  melted.  It  is  currently  assumed  that  snowflakes  have  an  ice  skeleton 
which  maintains  its  original  shape  during  most  of  the  melting  process  (Matsuo 
and  Sasyo,  1981a).  Fig.  11  shows  a  top  and  corresponding  side  view  of  a  snow¬ 
flake  melting  in  our  chamber.  This  ice  skeleton  can  easily  be  seen.  Matsuo  and 
Sasyo  postulate  that  this  ice  skeleton  collapses  to  a  water  drop  because  of 
surface  tension.  We  found  that  as  these  snowflakes  and  crystals  melted,  they 
frequently  disrupted.  Of  course,  caution  must  be  taken  in  interpreting  the 
disintegration  process  because  the  types  of  ice  crystals  we  could  catch  were 
limited  to  dendrites,  spatial  crystals  and  snowflakes  with  no  easily  identi¬ 
fiable  shape.  The  disintegration  of  ice  crystals  occurred  rapidly  in  the  final 
stages  of  melting  and  It  was  nearly  impossible  to  photograph  this  process  using 
the  present  camera  set  up.  However,  the  visual  observation  during  the  final 
stages  of  melting  clearly  supports  the  existence  of  the  disintegration  process. 

The  size  of  the  ice  crystal  affected  the  way  it  melted.  Crystals  less 


than  1  mm  in  diameter  were  easily  suspended  until  the  last  fractions  of  a  second 
of  their  crystal  form.  Then,  instead  of  falling  down  toward  the  bottom  of  the 
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SNOWFLAKE  DISRUPTED  1.5  s  LATER 

Fin  11.  Snowflake  reltinn  sentience. 


tunnel,  these  crystals  were  observed  to  take  a  violent  jump.  If  these  crystals 
had  transformed  into  spheres,  which  having  smaller  cross-sectional  areas,  they 
should  fall  faster.  In  our  observations,  this  downward  fall,  which  should  be 
a  uniform  increase  in  velocity,  was  the  exception.  These  small  crystals  jumped 
sideways  or  straight  up.  A  crystal  suspended  several  centimeters  from  the  tunnel 
wall  often  jumped  over  onto  it  and  appeared  instantly  as  a  water  droplet.  Static 
electricity  as  the  cause  of  this  was  ruled  out  because  of  our  test  with  the 
artificial  snow  crystal,  the  violent  nature  of  the  crystal  movement,  and  occur¬ 
ence  only  at  the  end  of  the  melting  process.  A  change  in  the  aerodynamic  con¬ 
dition  of  the  crystal  due  to  disruption,  could  cause  this  movement.  Total  dis¬ 
ruption  could  also  cause  the  pieces  of  the  crystal  to  be  swept  straight  up  and 
out  of  the  measurement  area. 

Disruption  also  offers  the  most  reasonable  explanation  for  the  snowflake 
behavior  in  the  final  seconds  of  its  melting  process.  The  snowflake  shown  in 
Fig.  11  was  easily  suspended  with  a  gradual  increase  in  the  tunnel  airflow 
throughout  the  pictures  taken.  However,  just  after  the  fourth  photograph,  it 
broke  up  and  fell  to  the  bottom  of  the  tunnel.  This  was  the  case  with  all  the 
snowflakes  we  suspended.  They  were  easily  suspended,  as  they  decreased  in  size 
by  melting,  up  to  a  certain  point  (twenty  seconds  being  the  longest  under  the 
temperature  and  relative  humidity  ranges  of  this  study),  whereupon  they  broke 
into  several  pieces,  with  the  larger  pieces  falling  to  the  bottom  of  the  tunnel 
and  the  smaller  ones  being  swept,  up  and  out  of  the  observation  section.  This 
occurred  as  if  there  were  a  small  explosion  in  the  snowflake.  The  larger  drops 
behaved  as  if  their  wings,  i.e.,  small  droplets  connected  by  the  ice  skeleton, 
were  lost,  falling  out  of  the  measuring  section  of  the  tunnel  before  the  operator 
could  react  to  increase  the  airflow  and  catch  them. 


Large,  flat,  symmetrical ,  dendritic  crystals  fell  flat  and  were  not  af¬ 
fected  by  this  disruption  phenomenon.  They  melted  slowly  from  the  outer  edges 
inwards  as  they  drifted  into  the  wall  (see  Fig.  12).  No  breakoff  of  a  portion 
of  the  crystal  was  observed.  This  action  of  drifting  into  the  wall  after  being 
stably  suspended  indicates  that  the  velocity  profile  suitable  for  a  snowflake 
or  a  crystal  is  not  sufficient  for  a  water  droplet. 

Environmental  Conditions 

Measurements  were  carried  out  under  the  following  combinations  of  relative 
humidity  (RH)  and  temperature  T. 

RH  (%) _ T  (°C) 

100  2 

3 

4 

90  3 

4 


3.2  Single  crystal  melting 

3.2.1  Results  of  measurement  and  discussion 

The  experimental  data  were  first  reduced  and  plotted  into  the  following 
4  relationships:  (i)  initial  diameter  vs  melting  time,  (il)  initial  area  vs 
melting  time,  (iii)  diameter  vs  time,  and  (iv)  velocity  vs  time. 

In  the  data  analysis,  a  clear  change  In  the  air  velocity  was  always  ob¬ 
served  on  the  chart  about  two  seconds  before  the  first  time  mark  of  photographing. 
Based  on  this  observation,  it  was  assumed  that  the  ice  crystals  were  being  sus¬ 
pended  for  two  seconds  before  the  first  photograph  was  taken.  We  estimated  the 
final  melting  time  on  the  recorder  chart  based  on  the  manner  the  stability  of 
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Fig.  12.  Photographs  of  a  dendritic  ice  crystal  melting  i 
the  vertical  wind  tunnel. 


the  ice  crystal  was  lost.  Also  the  data  for  veloctiy  and  diameter  changes  with 
time  were  classified  according  to  the  average  ice  crystal  size. 

Initial  Diameter  and  Area  vs  Melting  Time 

Figs.  13  and  14  show  the  observed  relationships  between  the  initial  diam¬ 
eter  dQ  and  melting  time  tm,  and  that  between  the  initial  horizontal  cross- 
section  area  AQ  and  tm-  The  longest  melting  time  occurred  for  non-dendritic 
ice  crystals  which  had  an  initial  diameter  of  about  3  mm  or  those  which  had 
an  initial  area  of  7.5  mm2.  In  the  figures,  the  non-dendritic  single  crystals 
showed  a  convex  upward  tendency  of  melting.  This  is  in  agreement  in  tendency 
with  the  treatment  for  ice  spheres  by  Mason  (1956)  and  Drake  and  Mason  (1966) 
and  that  which  we  have  done  for  non-dendritic  or  plate  crystals  as  we  shall 
discuss  later.  However,  they  are  different  from  dendrites  and  snowflakes,  whose 
melting  times  do  not  vary  appreciably  depending  on  their  sizes.  It  appears  that 
the  melting  of  dendrites  and  snowflakes  is  controlled  mostly  by  heat  conduction 
from  the  ventilating  air  through  their  structures.  It  also  seems  to  be  that 
dendritic  ice  crystals  with  initial  diameters  larger  than  3  mm  are  affected 
more  by  ventilation  than  those  smaller  than  3  mm.  The  AQ  vs  tm  plots  show  the 
same  tendency. 

Diameter  Change  With  Time 

Figs.  15  through  18  show  the  observed  diameter  d  vs  time  t  relationships. 

In  the  figures,  the  diameters  of  the  crystals  are  classified  into  seven  dif¬ 
ferent  groups;  they  are  1  (0  to  1  mm),  2  (1  to  2  mm),  3  (2  to  3  mm),  4  (3  to 
4  mm),  5  (4  to  5  mm),  6  (5  to  6  mm),  and  7  (6  to  7  mm).  It  is  natural  to  expect 
the  larger  the  ice  crystal,  the  longer  it  takes  to  melt.  As  can  be  seen  in 
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the  figures,  a  large  ice  crystal  does  show  a  quicker  change  in  diameter  than  a 
small  dense  ice  crystal.  Fig.  15  gives  slopes  only  for  the  groups  1,  2  and  3, 
but  it  suggests  that  the  groups  1  and  2  melt  slower  than  the  group  3.  Fig.  16 
shows  the  slopes  of  the  groups  2,  3  and  4  to  be  gentler  than  those  of  groups 
5  and  7.  This  seems  to  suggest  that  ice  crystals  in  the  0  to  3.0  mm  diameter 
range  melt  slower  than  those  in  the  3.0  to  7.0  mm  range.  Fig.  18  shows  clearly 
that  the  slope  of  the  group  4  is  steeper  than  either  the  group  2  or  3.  Whereas, 
Fig.  17  shows  an  exception  in  the  tendency  of  slope.  The  slope  of  the  group  2 
is  smaller  than  those  of  groups  3  and  4.  However,  these  data  indicate  generally 
that  slower  melting  ice  crystals  do  exist  in  small  diameter  ranges. 

Fall  Velocity  Change  With  Time 


Figs.  19  through  22  show  the  observed  fall  velocity  v  vs  time  t  relationship. 
In  these  figures,  it  may  be  seen  that  the  largest  total  change  in  fall  velocity 
occurs  in  general  for  smaller  crystals.  It  appears  that  large  dendrites  like 
the  group  7  not  only  fall  slowly  but  change  the  velocity  relatively  little  during 
the  melting. 

Figs.  21  and  22  show  the  variation  of  fall  velocity  with  time  for  ice  cry¬ 
stals  melting  under  relative  humidity  of  90%.  Although  the  relative  humidity  in 
these  cases  is  lower,  suggesting  slow-down  effect  of  melting  due  to  evaporation 
cooling,  the  initial  ice  crystal  conditions  are  different  compared  with  those 
in  Figs.  19  and  20.  Due  to  large  variation  of  crystal  conditions  involved  in 
the  measurements,  it  is  not  possible  to  clearly  isolate  the  effect  of  tempera¬ 
ture  or  relative  humidity  on  the  fall  velocity  variation  during  the  melting  from 
the  sample  data  analysis. 


Fig.  13.  Initial  diameter  of  single  Ice  crystals  d  vs  melting  time 
t_  for  all  the  environmental  conditions.  0  Details  of  the 

in 

environmental  effects  may  be  Interpreted  from  the  following 
figures. 


Initial  horizontal  cross-sectional  area  of  single  Ice  crystals 

vs  melting  time  tm  for  all  the  environmental  conditions, 
o  m 

Details  of  the  environmental  effects  may  be  Interpreted  from 
the  following  figures. 


Fig.  18.  Change  of  crystal  diameter  d  with  time  t  under  a 
relative  humidity  of  90%  and  a  temperature  of  4°C. 
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Fig.  20.  Change  of  crystal  fall  velocity  v  with  time  t  under  a 
relative  humidity  of  100%  and  a  temperature  of  3°C. 
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Fig,  21.  Change  of  cry-+al  fall  velocity  v  with  time  t  under 
relative  humidity  of  90%  and  a  temperature  of  3°C. 


Having  analyzed  the  data  for  single  ice  crystal  melting  in  a  vertical 
wind  tunnel  under  constant  temperature  and  relative  humidity,  we  found  that 
there  are  eight  groups  of  data  which  are  useful  to  construct  the  regression 
equations  for  the  time  evolution  of  ice  crystals.  The  conditions  of  the  tem¬ 
perature  and  relative  hunidity  of  each  group  are  listed  in  Table  1,  along  with 
the  initial  diameter. 


Table  1  Environmental  conditions  and  initial  Ice  crystal  diameter 
used  for  regression  analysis  of  single  crystal  melting. 


Group 

Temperature  (°C) 

Relative  humidity  (%) 

Initial  Diameter  (mm) 

A 

2 

100 

2.84 

B 

3 

100 

3.57 

C 

3 

100 

2.52 

0 

3 

100 

2.16 

E 

3 

90 

1.6 

F 

3 

90 

2.62 

G 

3 

90 

3.14 

H 

4 

90 

3.15 

The  regression  equations  constructed  for  each  group  of  data  are  as  follows 
for  the  diameter  (mm) 

d  «  dQ  +  dit  +  d2t2  +  d3t3  +  d4t4  +  d5t5,  (1) 

where  dQ  through  d&  are  constants. 

For  the  fall  velocity  (cm  s"*) 

v  -  vQ  +  Vjt,  (2) 

where  vQ  and  Vj  are  constants.  From  the  data  of  each  group,  we  have  determined 
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the  coefficients  In  (1)  and  (2)  used  to  draw  regression  lines  In  experimental 
data.  For  the  group  A,  the  regression  lines  for  the  diameter  and  the  fall 
velocity  as  a  function  of  time,  are  plotted  respectively  In  Figs.  15  and  19. 
For  groups  B,  C  and  0,  the  regression  lines  are  plotted  in  Figs.  16  and  20, 
respectively.  For  groups  E,  F  and  G,  they  are  plotted  in  Figs.  17  and  21 
and  for  group  H  in  Figs.  18  and  22. 

The  regression  equations  describing  the  change  of  diameter  and  fall 
velocity  of  the  snow  crystals  can  then  be  expressed  on  average  as 

d  -  d0  U  +  V  +  c2t2  +  c3t3  +  c4t4  +  V5)*  Ol 


V  =  dQ  (Cg  +  Cyt ) ,  (4 

where  dQ  is  the  initial  diameter.  The  constants  in  (3)  and  (4)  are  shown  In 
Tables  2  and  3. 


Table  2  Constants  for  the  regression  equation  describing  the 
diameter  of  snow  crystals  at  different  temperatures 
and  relative  humidities. 


T 

(°c) 

RH 

(X) 

Cj  (s-1) 

C2 

(s'2) 

C3  (s'3) 

c4  (s'4) 

c5  (s'5) 

2 

100 

-0.12 

5.27 

x  10’3 

0.91  x  10“6 

-0.33  x  10”4 

0.22  x  10" ! 

3 

100 

-0.07 

0.96 

x  10'2 

-0.24  x  10”2 

0.3  x  10”3 

-0.12  x  10' 

3 

90 

-0.06 

0.24 

x  10”2 

0.41  x  10’3 

-0.1  x  10”3 

0.56  x  10” 

4 

90 

-0.1 

0.07 

x  10'1 

0.3  x  10”2 

-0.035  x  10”  2 

0.15  x  10’ 

Table  3  Constants  for  the  regression  equation  describing  the 

fall  velocity  of  snovf  crystals  at  different  temperatures 
and  relative  humidities. 


Temperature  (°C) 

Relative  Humidity  (%) 

C6 

C7 

2 

100 

38.66 

1.17 

3 

100 

40.25 

1.27 

3 

90 

35.19 

1.45 

4 

90 

41.74 

0.6 

Furthermore,  we  can  generalize  the  equations  considering  the  environ¬ 


mental  temperature  =  Tq  +  AT,  Tq  being  the  melting  point  of  ice. 


d  =  dQ  (1  +  ajATt  +  o^ATt2  +  a^ATt3  +  a^ATt4  +  a^ATt^), 


(5) 


and 


v  3  do  <6o  +  61AT+  62ATt*' 


(6) 


where  3Q  through  e2  are  constants.  The  constants  in  (5)  and  (6)  are  given  in 
Tables  4  and  5.  We  can  notice  that  for  AT  *  0°C,  the  snow  crystal  will  keep 
its  size  with  constant  fall  velocity. 


Table  4  Constants  for  the  generalized  equation  describing  the  diameter 
of  snow  crystals  at  different  relative  tumidities. 


RH 


a 


(%)  (°c'Vl)  (°c"Vl) 


(°C’1s"J) 


rev*) 


(°c"V5) 


100  -0.35  x  10"1  0.3  x  10"2 


-0.56  x  10“3  0.64  x  10‘4 


•0.25  x  10 


-5 


90  -0.23 


-0.89  x  10* 


0.53  x  10"3  -0.69  x  10’4 


0.31  x  10* 
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Table  5  Constants  for  the  generalized  equation  describing  the  fall 
velocity  of  snow  crystals  at  different  relative  humidities. 


RH 

e2 

<*) 

(s'1) 

rc-v2) 

(•c-v*) 

100 

39.45 

0.06  x  10"1 

0.47 

90 

38.47 

0.13 

0.27 

3.3  Snowflake  melting 

3.3.1  Results  of  measurement  and  discussion 

The  previously  mentioned  problems  of  suspending  snowflakes  reduced 
the  number  of  cases  observed  to  42: 

Relative  Humidity  (%)  Temperature  (°C)  Number  of  Cases 


Snowflakes  were  Initially  identified  from  the  photographs  taken.  From  the 
observations,  particles  larger  than  2.5  and  1.0  mm  (a  and  c  dimensions) 
were  considered  snowflakes  and  Included  In  this  study. 

In  analyzing  the  data,  several  assumptions  were  made.  As  mentioned 
above,  an  enveloping  ellipsoid  was  usei?  to  represent  the  size  of  a  snow¬ 
flake.  No  attempt  was  made  to  compensate  for  the  amount  of  the  ellipsoid 
filled  by  the  outer  dimensions  of  the  snowflake.  Therefore,  the  snowflakes 
were  all  assumed  to  have  been  in  the  chamber  envi ronment  for  two  seconds  by 
the  time  the  first  photograph  was  taken  as  in  the  case  of  single  crystals. 
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Finally,  the  melting  time  had  to  be  estimated  due  to  the  disruption  phenome¬ 
non,  along  with  the  apparent  inability  of  our  air  velocity  profile  to  sus¬ 
pend  change  of  fall  velocities  for  like  groups  under  the  same  relative 
humidity  conditions  but  different  temperatures  (Figs.  26  and  27).  The  velo¬ 
city  increase  appears  to  be  faster  under  the  higher  temperature  (3°C)  because 
the  drop  is  melting  faster,  thus  its  cross-sectional  area  becomes  smaller 
more  rapidly. 

Fall  Velocity  Change 

Fig.  23  shows  the  fall  velocity  change  with  time  for  six  area  catagories 
of  snowflakes.  These  areas  represent  maximum  horizontal  cross-sectional 
areas,  perpendicular  to  the  direction  of  fall.  All  conditions  observed  are 
included  in  the  figure.  With  the  exception  of  the  15-20  and  25-30  mm*  area 
groups,  flakes  with  larger  areas  seemed  to  fall  faster  than  those  with 
smaller  areas.  A  comparison  of  the  fall  velocities  of  equal  area  groups 
under  the  same  temperature  but  different  relative  humidities  suggests  some 
interesting  information.  Fig.  24  shows  fall  velocities  under  90%  relative 
humidity  and  temperature  of  4°C.  Fig.  25  gives  fall  velocities  for  the 
same  temperature  but  with  a  relative  humidity  of  100%.  The  fall  velocity 
Increase  for  the  same  size  snowflake  appears  to  be  faster  in  the  100%  rela¬ 
tive  humidity  environment  than  in  the  90%.  This  may  be  explained  by  noting 
that  in  an  environment  of  90%  relative  humidity,  evaporative  cooling  or 
lowered  condensation  heating  reduces  the  rate  of  melting,  therefore  the  rate 
of  velocity  increase.  Whereas,  in  the  saturated  environment,  melting  occurs 
without  this  effect,  so  melting  proceeds  relatively  rapidly  causing  faster 
velocity  change.  This  tendency  is  seen  also  in  the  change  of  fall  velocities 
for  like  groups  under  the  same  relative  humidity  conditions  but  different 
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Fall  velocity  vs  time  at  4°C,  90%  RH. 


t  (s) 


Fall  velocity  vs  time  at  4°C,  100%  RH. 


temperatures  (Figs.  26  and  27).  The  velocity  Increase  appears  to  be  faster 
under  the  higher  temperature  (3°C)  because  the  drop  Is  melting  faster,  thus 
Its  cross-sectional  area  becomes  smaller  more  rapidly.  This  area  should  be 
studied  further. 

Shape  Change 

Fig.  28  shows  the  various  area  groups  to  be  decreasing  towards  final 
values  (the  areas  of  the  water  droplets).  Observations  proved  that  snow¬ 
flakes  melted  at  such  a  rate  that  their  suspension  was  controllable,  until 
they  disrupted.  This  disruption  size  was  not  predictable.  Because  this  dis¬ 
ruption  occurred  so  quickly,  a  photograph  could  not  be  obtained  for  analysis. 
For  analysis  of  the  data  shown  in  Figs.  28  thorugh  32,  one  must  be  cautious 
considering  the  disruption  phenomenon  that  exists  frequently  at  the  end  of 
the  measurement. 

Melting  Rates 

Examination  of  melting  time  versus  the  initial  diameter  of  the  snowflakes 
In  Fig.  33  yields  some  interesting  information.  Current  theory  of  snowflake 
melting  (Matsuo  and  Sasyo,  1981a)  predicts  a  curve  which  is  concave  upwards 
(by  treating  the  snowflake  as  a  solid  object),  while  our  data  shows  a  convex 
curve.  The  melting  time  increases  with  the  initial  size,  but  there  is  no 
appreciable  time  difference  between  sizes.  This  may  be  because  snowflakes 
are  of  a  porous  structure,  not  solid,  thus  allowing  the  air  to  pass  through 
them.  This  air  passing  through  the  structure  as  It  falls  makes  the  heat 
transfer  distance  from  the  passing  air  nearly  as  small  as  that  for  a  small 
Ice  crystal.  Thus  the  crystals  which  make  up  the  snowflake  can  melt  somewhat 
Independently  of  each  other.  The  larger  the  snowflake,  the  larger  the  Ice 
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Fig.  28.  Horizontal  cross- 


Fig.  29.  Horizontal  cross-sectional  area  A  vs  time  t  at  4°C,  90%  RH. 


Horizontal  cross-sectional  area  A  vs  time  t  at  2°C,  100%  RH 
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Fig.  31.  Horizontal  cross-sectiona 


Fig.  32.  Horizontal  cross-sectional  area  A  vs  time  t  at  4°C,  100%  RH. 


Fig,  33.  Melting  time 


skeleton,  which  would  make  this  process  more  Inefficient  thereby  giving  the 
slightly  longer  melting  time.  The  effect  of  relative  humidity  on  the  melting 
time  of  snowflakes  was  studied  by  Matsuo  and  Sasyo  (1982b  and  c).  Our  data 
shown  In  Figs.  34  and  35  are,  however,  not  sufficient  to  check  the  effect. 

3.3.2  Regression  analysis 

There  are  7  groups  of  snowflake  melting  data  which  are  useful  for  the 
regression  analysis  to  describe  the  time  evolution  of  snowflake  behaviors. 

The  conditions  of  the  temperature  and  relative  hunidity  of  each  group  are 
listed  in  Table  6,  along  with  the  intitial  area  and  initial  diameter  deter- 
minted  by  the  regression  analysis. 


Table  6  Environmental  and  initial  snowflake  conditions  used  for  re¬ 
gression  analysis  of  snowflake  melting. 


Group 

T 

<°C) 

RH 

(*) 

Ao 

(mmz) 

do 

(mm) 

A 

2 

100 

33.23 

6.51 

B 

2 

100 

22.38 

5.34 

C 

2 

100 

17.7 

4.75 

D 

3 

100 

17.51 

4.72 

E 

3 

100 

40.2 

7.15 

F 

4 

100 

14.19 

4.23 

G 

4 

100 

8.38 

3.16 

For  groups  A,  B  and  C,  the  regression  lines  of  the  change  of  area,  and 
fall  velocity  with  respect  to  time  are  plotted  in  Figs.  30  and  26,  respec¬ 
tively.  For  groups  0  and  E»  same  plots  are  shown  In  Figs.  31  and  27.  For 
groups  F  and  G,  they  are  shown  In  Figs.  32  and  27. 
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Assuming  that  the  surface  of  snowflake  is  spherical,  we  can  express  the 
cross-sectional  area  A  *  ^d2,  where  d  is  the  diameter.  The  regression  equa¬ 
tions  describing  the  time  variation  of  diameter  and  fall  velocity  of  the 
snowflakes,  therefore,  are  expressed  as: 

d  =  dQ  (1  +  Cjt  +  C2t2  +  C3t3  +  C/  +  C5t5),  (7) 

and 


v  -  d0  (v0  +  Vjt),  (8) 

where  dQ  is  the  initial  diameter.  The  constants  in  (7)  and  (8)  are  listed 
In  Tables  7  and  8. 

Table  7  Constants  for  the  regression  equation  describing  the 
diameter  of  snowflakes  at  different  temperature  and 
relative  humidity. 


T 

RH 

Cx 

c2 

c3 

c5 

(°C) 

(*) 

(s’1) 

(s’2) 

(s’3) 

(S-) 

(s’5) 

2 

100 

-0.13 

0.62  x  10’2 

0.25  x  10  "2  0.58  x  10"3 

-0.33  x  10’4 

3 

100 

-0.16 

0.16  x  10’1 

0.6  x  10"2 

-0.41  x  10"3 

0.29  x  10’5 

4 

100 

-0.12 

-0.55  x  10’2 

0.13  x  10’2  -0.64  x  10’4 

0.46  x  10’5 

Furthermore,  we  can  generalize  these  equations  considering  the  effect 
of  environmental  temperature: 

d  =  dQ  (1  +  ojATt  +  ct2ATt2  +  cijATt3  +  a4ATt4  +  c^ATt5),  (9) 


and 

v  -  dQ  (B0  +  SjAT  +  B2ATt), 


(10) 


71 


Table  8  Constants  for  the  regression  equation  describing  the  fall 
velocity  of  snowflakes  at  different  temperature  and 
relative  humidity. 


T 

(°C) 

RH 

{*) 

Vo 

(s-1) 

V1 

(S-8) 

2 

100 

16.91 

1.7 

3 

100 

18.64 

1.59 

4 

100 

29.23 

2.47 

to  describe  the  evolution  of  melting  snowflakes  under  constant  temperature 
AT  (°C)  and  100%  RH.  The  constants  in  (9)  and  (10)  are  given  in  Tables  9 
and  10.  It  is  clear  that  when  AT  *  0°C,  the  snowflake  will  keep  its  size 
with  constant  fall  velocity. 


Table  9  Constants  for  the  generalized  equation  describing  the 
diameter  of  melting  snowflakes. 


RH 

(*) 

ai 

rrv1) 

a2 

(°C”V8) 

°s 

(°C"VS) 

(°C‘V*) 

a5 

<°rv5) 

100  -0.042  -0.2  x  10"2  0.63  x  10'3  -0.11  x  10“4  -0.13  x  10’5 


Table  10  Constants  for  the  generalized  equation  describing  the 
fall  velocity  of  melting  snowflakes. 


RH  (%) 

0,(  s‘l) 

f^rc-v*) 

B2(°Cls”2) 

100 

21.59 

0.43 

0.62 
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3.3.3  Factor  e 


For  the  melting  of  snowflakes,  it  may  be  assumed  that  the  heat  transfer 
rate  to  the  surface  of  snowflake,  dQ/dt,  is  equal  to  the  sensible  heat  con¬ 
duction  from  the  environment  and  the  latent  heat  transfer  by  condensing 
water  vapor  on  the  surface.  For  a  spherical  snowflake, 

=  4irr(fhKAT  +  L^DAp),  (11) 

where 


r:  radius 

f^:  ventilation  coefficient  of  snowflakes  for  heat  diffusion 
K:  thermal  conductivity  of  air 

AT:  temperature  difference  between  the  environment  and  snowflake 
(AT  =  T.  -  Tw) 

L  :  latent  heat  of  condensation  of  water  vapor 

v» 

f  :  ventilation  coefficient  of  snowflakes  for  water  vapor  diffusion 
D:  coefficient  of  molecular  diffusion  of  water  vapor  in  air 

Ap:  difference  between  water  vapor  density  in  environment  and 
equilibrium  water  vapor  density  at  snowflake  surface. 

(Ap  "  pw,sat^ 

It  should  be  pointed  out  that  the  thermal  conductivity  of  water  is  more 
than  20  times  larger  than  that  of  the  air,  and  the  thickness  of  water  layer 
is  much  smaller  than  that  of  air  layer  through  which  heat  conduction  proceeds. 
Therefore,  we  assume  that  Tw  =  0°C  and  t  =  p0,  where  PQ  is  the  satura¬ 
tion  vapor  density  at  T  *  0°C.  Furthermore,  we  assume  that  for  a  non- 
spherical  snowflake,  (11)  is  appllable  with  a  factor  e  on  the  right-hand 
side  where  e  is  the  factor  comprising  all  the  effects  due  to  the  non-spherlcal 
shape. 

Also,  we  may  assume  that  the  total  amount  of  heat  transferred  from  the 
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air  stream  for  given  time  interval  At,  is  totally  consumed  by  melting  the  ice 
existing  In  the  surface  layer  of  Ar  thickness.  Therefore,  (11)  is  trans¬ 
formed  to 


-Lf4Trr2p^Ar  ■  4irre(fhKAT  +  LcfvDAp)  At,  (12) 

where 

L^:  latent  heat  of  fusion 
p.j :  density  of  ice  skeleton. 

As  At  0,  (12)  can  be  rewritten  as 

tfh“T  *  W*'-'-  (13> 

Therefore,  e  can  be  determined  from  (13),  with  the  ventilation  coefficients 
being  expressed  as 

f  -  fh  3  fv  *  1  +  °’ZA  (14) 

(Matsuo  and  Sasyo,  1981),  where  Re  Is  the  Reynolds  number. 

Obviously,  (14)  Indicates  that  the  value  of  e  changes  with  time.  Hence, 
applying  the  regression  equations  about  diameter  and  fall  velocity  described 
earlier,  we  may  obtain 

e  *  dQ  (e0  +  Ejt1  +  e2t2  +  Cjt3  +  e4t4  +  e5t5).  (15) 

The  constants  in  (15)  fbr  snowflakes  are  listed  in  Table  11. 

Further,  we  may  generalize  the  equation 

e  -  dQ  AT  (d0  ♦  bjt  +  b2t2  +  b3t3  +  b4t4  +  bgt5),  (16) 

for  describing  the  change  of  e  of  the  snowflake  melting  with  constant  tem¬ 
perature  AT  (°C)  and  100*  Relative  Humidity.  The  constants  In  (16)  are  given 
in  Table  12. 
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Table  11  Constants  in  equation  {15)  for  describing  the  e  of  snowflakes. 


T 

RH 

eo 

£i 

£3 

£5 

’C) 

(*) 

(mm-1) 

(mm"1s"1) 

(mm'V2) 

(mm'V3) 

(mm'V4) 

(mm'V*) 

2 

100 

4.09 

-0.45 

-0.19  xlO'1  -0.16  xlO'2 

-0.59x10' 3 

0.99  xlO"4 

3 

100 

3.45 

-0.16 

-0.9 

0.11  x 10" 1 

0.34  xlO'2 

-0.3  xl0‘* 

4 

100 

1.52 

-0.17 

0.7  x  10"2  -0.61  x  10"3 

-0.65  xlO'3 

0.74  x JO'4 

Table  12  Constants  in  equation  (16)  for  the  generalized  description  of  t. 


RH 

(%) 

bo 

(n*n°C)'1 

bx 

(mm°Cs)"1 

b2 

(mm°Cs2)'1 

lmrn°Cs3)”1 

K 

(mm°Cs4Ti,:i 

bS 

(mm°Cs5)“1 

100 

0.85 

-0.71  xlO'1 

-0.2  xlO'1 

0.99  x 10" 3 

0.22  xlO'3 

-0.14  xlO'4 

4.  Theories  of  Ice  Crystal  and  Snowflake  Melting 

Mhen  ice  melts,  the  latent  heat  of  fusion  is  absorbed  by  the  ice  and 
the  heat  Is  normally  supplied  from  outside.  So,  melting  proceeds  inwards 
from  the  outer  edges  and  surfaces  of  the  Ice. 

In  melting  of  snow  crystals  and  snow  flakes  during  their  fall  through 
zones  of  temperature  higher  than  0°C,  the  manner  the  heat  is  supplied  to 
the  ice  surfaces  controls  the  melting.  Heat  may  be  transported  to  the  ice 
surfaces  through  the  air  normally  in  the  following  three  mechanisms;  radiative 


transfer,  sensible  heat  conduction  and  latent  heat  transfer  by  water  vapor 
diffusion.  While  the  atmospheric  ice  melting  proceeds,  the  cloud  above  and 
the  falling  melted  crystals  below  the  melting  layer  sustain  a  radiative 
balance  in  the  zone.  Under  such  a  balance  at  a  temperature  as  low  as  0°C, 
the  radiative  heat  transfer  to  the  melting  layer  may  be  ignored.  The  heat 
transfer  to  the  surface  of  melting  ice  in  the  atmosphere  may  therefore  be 
described  by  two  other  mechanisms,  i.e.,  sensible  and  latent  heat  transfers. 

If  an  ice  crystal  is  placed  in  quiescent  warm  air  without  fall,  a  tem¬ 
perature  field  will  establish  around  the  crystal  and  sensible  heat  conduc¬ 
tion  occurs  from  the  environment  to  the  crystal  through  the  field.  The 
latent  heat  transfer  takes  place  in  an  identical  manner  but  the  direction 
depends  on  the  condition  of  the  environment.  If  the  environment  is  dry, 
heat  removal  takes  place  due  to  evaporation  cooling.  These  processes  are 
reasonably  well  understood  with  respect  to  ice  crystal  growth  under  super¬ 
cooled  environment  although  their  directions  may  be  partially  or  totally 
reversed. 

Ice  melting  in  the  atmosphere  takes  place  while  ice  crystals  fall. 

The  fall  results  In  sheared  air  motions  (convection)  around  the  crystals, 
thereby  deforming  the  temperature  and  water  vapor  fields  of  quiescent  air. 
The  fall  motion  of  an  Ice  crystal  modifies  the  field  most  extensively  at  the 
distal  points  If  the  crystal  shape  were  plate-like  or  prism-like.  The  dis¬ 
tal  points  are  for  this  reason  the  zones  of  fastest  Ice  crystal  melting. 

As  ice  melting  proceeds,  formed  water  should  show  a  tendency  of  taking  a 
shape  of  smaller  surface  area.  Since  the  shape  that  takes  the  smallest 
surface  area  for  a  given  volume  Is  a  sphere,  the  Ice  crystal  shrinks  from 
the  distal  points.  This  shape  change  of  Ice  crystal  during  the  melting 
leads  to  an  Increase  In  the  fall  speed.  The  fall  speed  Increase  In  turn 


alters  the  convective  air  motion  around  the  crystal  and  further  modifies 
the  accompanied  temperature  and  vapor  fields. 

The  fall  speed  change  during  the  melting,  although  It  is  beyond  the 
scope  of  present  study,  will  provide  Important  leads  to  estimation  of  ice 
and  water  content  variation  as  a  function  of  altitude  and  to  interpretation 
of  radar  reflectivity  data  in  the  bright  band. 

We  shall  next  attempt  formulation  of  this  atmospheric  ice  melting 
process  first  with  single  crystals  and  then  with  snowflakes. 

4.1  Single  ice  crystal  melting 

4.1.1  Description  of  the  problem  and  the  approach 

There  are  two  basic  shapes  of  atmospheric  single  ice  crystals,  i.e., 
plate  and  prism.  In  order  to  permit  a  smooth  theoretical  treatment  of 
various  factors  in  atmospheric  ice  melting,  it  is  necessary  to  simulate 
the  two  basic  ice  crystal  shapes  with  simple  geometrical  bodies.  A  common 
practice  in  this  regard  is  to  use  rotational  spheroids,  oblate  for  plate 
crystal  and  prolate  for  prism  crystal.  For  simplicity,  from  now  on  we 
shall  discuss  oblate  spheroid  only.  Treatment  with  a  prolate  spheroid  can 
be  easily  obtained  in  an  identical  manner  as  we  shall  discuss  below. 

First  of  all,  it  is  possible  to  treat  the  problems  of  heat  conduction 
and  vapor  diffusion  to  and  from  the  ice  crystal  in  quiescent  air  environ¬ 
ment,  and  then  to  add  the  effect  of  convective  motion  of  air  around  the  cry¬ 
stal  which  is  caused  by  the  fall.  In  stationary  air,  the  heat  conduction 
and  vapor  diffusion  around  the  rotational  oblate  spheroid  can  be  conveniently 
expressed  by  a  spheroidal  coordinate  as  shown  in  Fig.  36.  Since  the  effects 
of  condensation  and  thermal  accommodation  coefficients  are  negligible  for  the 
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ISOBAR  OR 
ISOTHERM 


Fig.  36.  Rotational  oblate  coordinate  for  heat  conduction  and 

vapor  diffusion  around  a  stationary  Ice  crystal  In  air. 

Ice  crystal  size  In  question,  the  Isotherms  of  heat  conduction  and  the  Iso¬ 
bars  of  vapor  diffusion  correspond  to  the  set  of  confocal  spheroids  outside 
the  Ice  crystal.  Similarly,  the  flow  lines  of  heat  and  vapor  are  Included 
In  the  set  of  confocal  hyperboloids  of  two  sheets.  Since  the  coordinate  Is 
an  orthogonal  system,  the  set  of  spheroids  and  the  set  of  hyperboloids  cross 
each  other  at  right  angle.  An  analytic  solution  for  heat  or  water  vapor  flux 
at  an  arbitrary  point  of  this  crystal  surface  Is  available. 

The  thermal  flux  at  the  distal  points  of  Ice  spheroid  Is  the  largest  and 
as  the  crystal  melts,  the  surface  tension  forces  the  formed  water  to  gather 
mostly  at  proximal  areas,  and  let  the  water  lenses  sandwich  the  crystal  (see 
Fig.  37).  Under  such  a  double  system,  the  heat  transfer  by  the  sensible  and 
latent  heat  fluxes  suddenly  becomes  complex.  There  exists  a  clear  need  of 
proper  simplification  for  the  process.  By  examining  the  thermal  conductivity 
In  air  (Ka),  that  In  Ice  (K^)  and  that  In  water  (1^); 


<T„‘  pj 


WATER. 


ICE  CRYSTAL 
(0°C) 


4Jw*  pw,sat^ 


WATER 


Fig.  37.  A  double  rotational  oblate  spheroid  model  for  partially  melted 
ice  crystal.  Both  water  and  ice  spheroids  are  in  contact  at 
the  distal  point  due  to  the  surface  tension  effect  of  water. 

Ka  =  5.73  x  10'5cal  (cm  *  s  *  °KP  (0°C) 

Ki  =  4.99  x  10"3 

-  1.40  x  10' 3  "  ’*  (20°C) 

1.52  x  10' 3  "  H  (50°C) 


it  may  be  seen  that  1^  >  20  Ka  at  0°C.  This  is  to  say  that  the  heat  conduc¬ 
tion  is  slower  in  air  under  steady  state  than  in  water.  So,  in  addition  to 
simulating  water  lenses  sandwiching  melted  ice  with  a  rotational  oblate 
spheroid  (see  Fig.  37),  we  assume  that  this  double  spheroid  consisting  of 
water  and  ice  spheroids  behaves  In  the  same  manner  as  an  identical  Ice 
spheroid  as  far  as  heat  and  water  vapor  fields  of  the  surrounding  air  are 
concerned . 
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When  this  crystal  falls,  the  convective  air  motion  further  changes  the 


heat  conduction  to  the  spheroid.  However,  the  overall  effect  of  this  air 
ventilation  is  known.  It  is  clear  that  the  ventilation  effect  is  zero  at  the 
proximal  points  and  is  largest  at  the  distal  points.  For  this  reason,  we  pro 
ceed  to  assune  that  the  ventilation  effect  is  proportional  to  the  heat  flux 
at  the  distal  points. 

With  these  assumptions,  we  shall  attempt  to  estimate  the  rates  of 
changes  of  major  semiaxis  a,  and  minor  semiaxis  of  water  spheroid  c  and  that 
of  ice  spheroid  c'  (see  Fig.  37). 


4.1.2  Melting  of  plate-like  ice  crystal 
(a)  Formulation  of  the  theory 


When  an  Ice  crystal  melts,  the  rate  of  heat  absorption  relates  to  the 
rate  of  melting  in  the  following  manner 


dQ 

3F 


W  <at>- 


(17) 


where  Q  is  the  heat,  t  the  time.  If  the  latent  heat  of  fusion  and  m  the  mass 
of  melted  ice.  Since  the  heat  is  supplied  by  two  mechanisms,  l.e.,  sensible 
and  latent  heat  (vapor  diffusion)  transfers, 

aT  "  *dQs  +  dQl*’ 

where  subscript  s  stands  for  sensible  and  1  for  latent.  Using  analogy  to  the 
electrostatic  field  around  a  conductor  condensor  of  oblate  shape,  (analogy 
was  also  used  for  the  temperature  and  vapor  fields  of  a  growing  Ice  crystal 
In  a  quiescent  air  environment),  we  have 

a?  * « [  K.  <T.  -  V  *  DLc  -  <w>  ]  •  (19> 
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where  C  is  the  electrostic  capacitance  of  the  outer  water  oblate  spheroid 
enveloping  the  inner  oblate  ice  spheroid,  being  given  as 
ae 

c  ■  •  (20) 

where 

e  -  (1  -  c2/a2)*. 


Since  T^  and  pw  sa^  are  not  known,  (19)  is,  strictly  speaking,  not  a  solu¬ 
tion.  However,  as  we  have  examined  above,  the  thermal  conductivity  of  water 
is  more  than  20  times  larger  than  that  of  air.  Furthermore,  the  thickness 
of  the  water  layer  is  much  smaller  than  that  of  the  air  layer  through  which 
heat  conduction  must  proceed.  So,  it  may  be  reasonable  to  assume  that 

Tw  -  V  and  pw,sat  ~  V 

where  p0  is  the  saturation  vapor  density  at  T  ■  T0  and  TQ  is  the  melting 
point  of  ice  or  0°C.  Since  ice  melting  normally  takes  place  in  air  whose 
temperature  is  not  far  above  0°C  and  whose  relative  humidity  is  high,  and 
air  flow  condition  around  the  falling  and  melting  ice  crystals  is  gentle, 
this  assumption  is  not  expected  to  create  a  large  amount  of  error. 

So  far,  we  have  obtained  the  expression  for  heat  flux  which  is  used 
for  ice  melting  when  the  ice-containing  oblate  water  spheroid  is  placed  in 
quiescent  air.  When  the  spheroid  falls  through  the  air,  the  air  motion 
changes  the  thermal  and  vapor  fields  around  it  (see  Fig.  38).  This  ventila¬ 
tion  affects  the  fields  differently  but  there  are  reasons  to  believe  the 
difference  is  small  and  the  normal  practice  is  to  set  both  corrections  same. 
Let  us  describe  the  ventilation  correction  by  f.  Then,  we  can  modify  (19) 
to  obtain 
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Fig.  38.  Numerically  computed  stream  function  arc!  vortlcity 
distribution  around  an  oblate  speroid  of  c/a  *  0.05 
at  (a)  Re  *  1.0  and  (b)  Re  *  20  (Pruppacher  and 
Klett,  1978). 


dQ 

at 


4wC  f  [  Ka  (T.  -  T0)  ♦  0LC  (p_ 


(22) 


The  above  equation  already  includes  Tw  -  TQ,  px  =  pQ  approximation,  and  con¬ 
sists  of  known  quantities  only. 

The  above  treatment  has  clarified  heat  transfer  to  a  ventilated  water 
spheroid  containing  a  melting  ice  spheroid.  Now,  we  try  to  understand  the 
shape  change  inside  the  water  spheroid  during  the  melting.  The  righthand 
term  of  (17)  may  be  expressed  as 


.  dm 
Lf  dt  ' 


,  dVL 

Lf°L3r 


(23) 


where  V,  and  p.  are  respectively  the  volume  and  the  density  of  liquid  water. 

L 

If  we  can  ignore  evaporation  loss  of  water,  the  volume  of  the  outer  water 
spheroid 


V 


+  Vc 


(24) 


where  V$  is  the  volume  of  the  ice  spheroid,  and 
4  9 

V  *  ^  Tra2c  , 


(25) 


4  2  - 

Vs  =  5  ™2c  , 

(see  Fig.  37).  Then  from  eqs.  (24),  (25)  and  (26) 

VL  =  V  -  Vs  =  ^ira2(c  -  c")  . 

Then 

dVL  *  dV  -  dV$  *  ^ir[2a(c  -  c'  )  da  +  a2(dc  -  dc'  )]  . 


(26) 


(27) 


(28) 


83 


From  mass  conservation  Ignoring  the  mass  loss  by  evaporation, 

Plvl  *  psvs  ■  V  {29) 

where  is  the  mass  of  ice  crystal  before  melting. 

It  is  known  that  the  sensible  heat  flux  to  the  surface  of  a  stationary 
spheroid  has  the  following  relationship 

(Heat  flux  on  the  a-axis)  a 

■  ■  —  --  ■  ■  — - -  ■  a  — 

(Heat  flux  on  the  c-axis)  c 

An  identical  relationship  holds  for  the  latent  heat  transfer  by  vapor  dif¬ 
fusion.  Therefore  the  above  relationship  is  valid  for  the  net  heat  flux. 
Since  these  heat  fluxes  are  used  to  melt  ice  crystal  at  corresponding  posi¬ 
tions. 


(Net  heat  flux  on  the  a-axis)  _  a  da  /,nx 

(Net  heat  flux  on  the  c-axls)  ”  c  3c"'  *  ' 

Note  that  although  da  corresponds  to  the  change  on  the  a-axis,  dc%  not  dc, 
does  that  on  the  c-axis  because  c'  represents  a  dimension  of  the  ice  crystal 
spheroid. 

With  the  convective  air  motion  around  the  falling  ice  crystal,  the  ven¬ 
tilation  effect  enhances  the  melting.  As  discussed  above,  the  ventilation 
effect  is  assuned  to  be  same  for  both  sensible  and  latent  heat  (vapor  dif¬ 
fusion)  conductions.  The  factor  f  relates  the  total  effect  of  ventilation 
to  the  stationary  heat  and  vapor  fields.  It  is  clear  that  there  is  no  ven¬ 
tilation  effect  on  the  c-axis  but  the  effect  is  at  Its  maximum  on  the  a-axis. 
The  effect  on  the  a-axis  must  be  proportional  to  (f  -  1).  Taking  a  propor¬ 
tionality  constant  Aj,  (30)  may  be  modified  to 
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af'  ‘  I  [  Ai  (f  - 1)  ♦  i  ]  •  on 

In  this  treatment,  we  leave  Aj  determination  for  experiments. 

We  now  have  three  relationships  for  variation  of  a,  c,  and  c',  l.e., 
eqs.  (28),  (29)  and  (30).  Eq.  (22)  adds  the  heat  flux  effect  or  Introduces 
time  variation  to  the  quantities  Involved.  So,  we  shall  solve  for  da,  dc 
and  dc'  using  eqs.  (28),  (29)  and  (30),  applying  a  fixed  time  interval  At 
to  Incorporate  the  net  heat  flux  effect.  In  order  to  solve  these  equations, 
we  convert  (29)  into  the  following  form, 

pLt^L  +  PsdVS  = 

or  inserting  (28)  and  the  differential  form  of  (10),  we  have 

2a  [  -  (pL  -  p$)c'  +  cpl  ]  da  +  a2pLdc  -  a2(pL  -  Ps)dc'  =  0.  (32) 

The  dV^  term  in  (28)  is  caused  by  the  net  heat  flux  in  (22).  So 

dV.  -  -ip  dQ,  (33) 

L  Lf 

and  from  (22) 

dQ  =  A2Cf,  (34) 

where 

*2  *  4"|>a<T»  '  V  +  °Lc<0.-po>]  «•  <35> 

In  other  words,  dQ  Is  the  amount  of  heat  used  to  melt  the  ice  crystal  in  At 
time.  This  introduction  of  time  interval  At  permits  numerical  computation 
as  a  function  of  time  allowing  Incorporation  of  effects  of  time  dependent 
variables.  Eq.  (28)  becomes 
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A2Cf  4 
\TfXr 


[  2a{c  -  c')da  +  a2dc  -  a2dc"J  . 


Using  eqs.  (31),  (32)  and  (36)  we  can  solve  for  da*  dc  and  dc"; 


A„A,C 


where 


3  "  7c 


~  At(f  -  1)  +  1  -  1 


dc  *  [’(pL  -  pS)c'  +  cpl]  [Al(f  “  1)  +  l]  “  (Pl-Ps*]* 


o  A«A«%C 


(b)  Computation  scheme 


With  these  equations,  it  is  now  possible  to  compute  melting  of  an  ice 
crystal  with  the  shape  of  rotational  oblate  spheroid  falling  in  air  with 
vertical  temperature  and  relative  humidity  distributions. 


COMPUTATION  SCHEME 


INITIAL  CONDITIONS 


a  -  a0,  c  »  cQ,  c1  -  c0,  T  -  T0,  p ,  p„.  v,  mQ 


:*N£ 


ASSOCIATED  CONDITIONS 


dQ  (34) 


da  (37),  dc  (39),  dc'  (40) 


a(=  aQ  +  da),  c(=  cQ+dc),  c'(=  cq  +  dc') 


V  +  dVL,  V$  +  dV$,  V  + 


n(a  +  da)(c+dc),  ir(a  +  da)(c'  +  dc'),  ir(a+da): 


V  +  dV 


4.1.3  Ventilation  coefficient  of  vapor  diffusion  and  heat  conduction 


According  to  Pruppacher  and  Klett  (1978),  the  simplified  ventilation 
coefficient  for  oblate  crystal,  f,  is  given  as 


f  »  1  +  0.142X2  +  0.054X,,ln(0.893X2), 


where 


X  ■  ^Sc^*  -  °'71  or  Re  -  0,63  *ith  NSc  "  °'71' 

where  Ngc  is  the  Schmidt  number  given  as 

N$c  *  ^  *  (v;  kinematic  viscosity), 

and  Re  is  the  Reynolds  number  given  as 


and 

f  *  0.937  +  0.178X,  (42) 

for  X  >  0.71  or  Re  >  0.63. 

When  X  il,  f  varies  linearly  with  X,  and  in  our  present  model,  X  is 
expected  to  fall  in  this  range.  These  authors  state  that  for  Re  <  100, 

f  *  0.86  +  0.28X  (X  i  1.0),  (43) 


agrees  with  experimental  data  to  within  8%.  So,  (43)  appears  to  be  our  best 
choice  for  the  ventilation  coefficient. 

For  the  terminal  velocities  of  oblate  spheroid  ice  crystals,  Pruppacher 
and  Klett  (1978)  list  the  following  empirical  relationship;  for  c/a  »  0.2, 

log  Re  *  Bq  +  BjX  +  B2X2,  (44) 

where  X  ■  log  (C0  Re2)  and  BQ  ■  1.3300,  Bj  •  1.0217,  B?  ■  -0.049018,  and  CD 
Is  the  hydrodynamic  drag  force  coefficient,  end  for  c/e  *  0.5, 

log  Re  «  B0  +  BjX  +  B^2  +  B3X3,  (45) 

where  BQ  -  -1.3247,  Bj  «  1.J396,  B2  ■  -0.047556,  B3  -  -0.002327.  For  an  ob¬ 
late  spheroid, 
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«D  Re! 


32a!c  (ps  -  pa)p  9  a«0psg 

S - - —  •  - - -  -  Mt  -  -  * 


where  pa  Is  the  density  of  air,  g  the  gravity  acceleration  and  n  the  dynamic 
viscosity  of  air.  They  Indicate  that  these  equations  are  applicable  for 
Re  <  100,  and  beyond  this  level  of  the  Reynolds  number  oscilliations  will 
occur  and  eventually  becomes  a  gliding- tumbling  motion. 


4.2  Snowflake  meltim 


Snowflake  melting  problem  has  been  recently  tackled  by  Matso  and  Sasyo 
(1981a).  Considering  uncertainties  involved  In  the  process,  their  theory 
appears  to  be  the  best  one  we  can  choose  at  this  moment.  So,  we  briefly 
review  it  below. 

The  rate  of  heat  transfer  by  the  sensible  heat  conduction  and  vapor 
diffusion  (latent  heat  transfer)  can  be  expressed  as,  similar  to  (22) 

3T  *  W**'7-  -  V  '  0Lc  •  W>]  •  (47> 


where  r  is  the  radius  of  snowflake,  e  a  factor  involving  the  effects  of 
various  physical  properties  of  the  snowflake  on  the  heat  transfer  rate, 
e  is  assumed  to  be  constant.  For  the  Reynolds  number  between  10  and  1800 

f  =  1  +  0.275  Npr3  Re*.  (48 

where  N-  is  the  Prandtl  number  being  given  as 
Pr 

\r  ■  ^a  •  (4! 

Since  Npr1  does  not  vary  much  In  our  range  of  atmospheric  conditions,  it 
is  taken  to  be  0.87.  This  reduces  (48)  to 


f  -  1  +  0.24  Re8  . 


According  to  Matsuo  and  Sasho,  during  the  melting  of  a  snowflake,  the  water 
penetrates  Inwards  leaving  the  surface  uncovered  with  water.  The  flux  of  the 
sensible  heat  and  that  of  the  latent  heat  match  the  rate  of  heat  absorption  by 
ice  melting,  or 


-  Lf4wr2pjdr  *  4irref(KaAT  +  DLcAp)At,  | 

where  is  the  density  of  ice-skeleton.  Eq.  (51)  with  simplifications 
discussed  above  reads 


-&=^Ki+o-24(IH- 


They  report  the  solution  of  (52)  as 

t  1  b(r  3/2  _  rs/2X 

1  a*RZTF  T  lro  r  ' 

-  W*(r#  -  r)  ♦  WrJ  -  r“>  -  )  .  (H) 

where  a  *  0.24  (2v/v)^  and  rQ  is  the  Initial  snowflake  radius.  Their  inte¬ 
gration  was  carried  out  under  the  condition  v  =  const.,  but  in  our  case,  v 
is  coupled  with  r  and  this  consideration  is  necessary  when  Integration  of 
(52)  is  performed. 

_  2 

In  the  range  10'  <  Re  <  300,  the  following  empirical  expression  may 
be  used  to  relate  the  fall  velocity  to  the  size  (Pruppacher  and  Klett,  1978) 

Y  -  BQ  *  BjX  +  ...  BgX6  ,  (54) 

where  8Q  ■  -0.318657  X  10\  B1  «  +0.992696,  B2  -  -0.153193  X  10"2,  B3  - 
-0.987059  X  10"\  B4  «  -0.578878  X  10‘»,  Bg  -  +0.855176  X  10“\  Bg  •  -0.327813 
X  10’ 5  and  X  »  in  (C^  Re2)  and  In  Re  *  Y. 


—  -  1 

5.  Mechanism  of  Bright  Band  Formation 

Our  vertical  wind  tunnel  studies  have  produced  some  Information  about 
variation  of  microphysical  parameters  during  the  melting  and  they  are  now 
represented  by  the  regression  equation  as  discussed  above.  Although  our 
knowledge  is  still  limited.  It  may  now  be  worthwhile  to  apply  it  to  a  real¬ 
istic  atmospheric  condition  to  see  the  contribution  of  Isolated  melting  particles 
to  the  radar  reflectivity.  Taking  this  as  a  stepping  stone  and  comparing  It 
with  observed  bright  band  data,  we  can  then  estimate  the  degree  of  contri¬ 
bution  from  other  factors.  Remembering  the  uncertainties  and  limitation  asso¬ 
ciated  with  our  data  and  present  radar  meteorology  In  general,  we  shall  proceed 
with  modeling  of  melting  layer  by  making  suitable  approximations. 

5.1  Model  of  melting  snowflake 

In  order  to  estimate  the  backscattering  from  a  melting  snowflake,  Ryde 
(1946)  considered  that  during  the  melting  process,  it  consisted  of  a  homo¬ 
geneous  mixture  of  Ice  and  water.  Austin  and  Bemis  (1950)  applied  the  method 
to  melting  snowflakes  to  compute  the  radar  reflectivity  in  the  melting  zone. 

However,  the  treatment  is  considered  to  be  questionable  (Battan,  1973). 

Clearly  the  effect  of  water  skin  on  the  melting  ice  has  to  be  considered. 

Whereas,  considering  the  wavelengths  of  microwave  commonly  used  In  radar  and 
the  thickness  of  the  skeletal  structure,  the  outer  dimension  of  a  melting 
snowflake  cannot  be  equivalent  to  that  of  solid  Ice  under  melting. 

The  effect  of  water  skin  existing  on  melting  spherical  Ice  particles 
was  estimated  by  Langleben  and  Gunn  (Battan,  1973)  based  on  the  work  by  Aden 
and  Kerker  (1951).  However,  the  effect  on  melting  snowflakes  Is  not  avail¬ 
able  and  Its  exact  computation  Is  beyond  the  scope  of  the  present  study. 
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There  are  uncertainty  and  variability  In  microphysical  and  environ¬ 
mental  conditions  as  well  as  limitation  on  our  knowledge  about  radar  back- 
scattering  from  melting  hydrometeors.  Considering  this  position,  we  have 
introduced  the  following  simplified  melting  snowflake  model  (see  Fig.  39). 


Fig.  39.  Melting  snowflake  model  involving  melting  ice  particles 
whose  oositions  are  fixed  with  the  space  coordinate. 


The  model  assunes  a  spherical  shape.  Initially,  it  has  a  radius  r^  and  an 
apparent  density  pf,  and  consists  of  nQ  number  of  ice  spheres  with  radius 
rpi.Q.  The  ice  spheres  are  fixed  with  the  space  coordinate  so  that  they 
cannot  change  the  positions.  As  the  melting  proceeds,  the  snowflake  melts 
from  outside  as  well  as  inside  and  allows  the  outer  ice  particles  to  disappear. 
The  water  formed  by  melting  is  assuned  to  distribute  equally  in  the  skin  of 
remaining  inner  ice  particles.  No  mass  exchange  is  considered  between  the 
snowflake  and  the  environment.  This  assumption  should  be  reasonable  in 
view  of  the  fact  that  melting  takes  place  normally  in  a  nearly  saturated 
environment  close  to  0°C  isotherm. 

The  ice  particles  are  assumed  to  change  the  radii  according  to  the 


relation 


where 


rp1  =  rp1oF» 


r“  rw 
rf  ■  rw  ' 


and  r  the  snowflake  radius  and  rw  the  radius  of  formed  water  drop  at  the 
end  of  melting. 

From  the  mass  conservation  assumption,  the  mass  of  an  ice  particle  with 


water  skin 


"p  ■  T,rfpf/n 

*  1”rpiPS  *  1’(rpi  -  rp1)pL’ 


where  n  is  the  number  of  remaining  ice  particles  in  the  snowflake  and  rpi 
and  r  the  radius  of  the  ice  particle  and  that  of  water-covered  ice  par- 

pw 

tide,  respectively.  Eq.  (57)  reduces  to 


nr*  - 

pw  PL  f 


S  ~  L'  r  3 
PL  Pi  * 


Snowflake  melting  from  outside  with  ice  particle  position  fixed  leads 


n/n0  *  3’irV|rrr^  =  r*/*^. 


Whereas,  the  initial  condition  is 


noKi°ps  *  l"rfpf 


r  s  ,  I  r* 

P10  VS  f  * 


Y>*'  A 


-|*cd 


Combination  of  eqs.  (55),  (58),  (59)  and  (60)  yields 


Pf(ps  ~  P[_) 
PLPS 


r*F(r) 


(61) 


To  incorporate  the  effect  of  the  water  skin  on  melting  ice  particles, 
we  assume 


o  =  o,  -  (cr.  -  ojG(r), 


(62) 


where  a  is  the  backscattering  cross  section  of  a  single  particle  and 


G(r) 


rJ  _  rJ 

ITT?' 

f  w 


(63) 


where  j  is  the  number  to  be  determined  from  the  computation  made  by  Langleben 
and  Gunn.  The  estimation  by  Langleben  and  Gunn  was  for  ice  particle  radius 
of  0.6  mm  and  the  wavelength  of  microwave  of  3  cm.  It  is  rather  insensitive 
to  the  wavelength  range  of  commonly  used  radar.  This  permits  us  to  apply 
their  results  to  particle  radii  around  0.6  mm.  For  constant  outer  radius  of 
water- coated  ice  particle. 


where 


1K|  =  . 


|k|*  •  |£rii  • 


(64) 


m  being  the  complex  index  of  refraction.  This  leads  to  the  relationship 

|K|!  -  |K,|’  -  (|K,|!  -  |g!)  G(r).  (65) 


Then,  since  the  radar  reflectivity 
n«  i*l2  V2  N, 

where  V  Is  the  volume  and  N  the  number  concentration  of  snowflakes,  from 
eqs.  (61),  (65)  and  (66) 
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i 


I 


H 


1  «  [l K, | ’  -  (|K,|!  -  |g»)  S(r)][r»  -  ^  r»F(r)]|,  (67) 

where  v  is  the  fall  velocity  of  snowflake. 

5.2  Numerical  simulation  of  the  melting  of  snowflakes  in  a  realistic 

environment 

Assuming  the  existence  of  a  realistic  atmospheric  environment,  we 
may  numerically  simulate  the  melting  of  either  snowflakes  or  ice  crystals 
in  it,  as  they  fall  through  the  region  where  the  air  temperature  is  higher 
than  ice.  Under  the  process,  the  change  of  the  size  of  snowflakes  or  ice 
crystal,  and  the  change  of  its  fall  velocity  may  be  calculated  according  to 
the  generalized  expressions  derived  from  data  obtained  in  the  vertical  wind 
tunnel  under  constant  temperatures.  We  shall  outline  the  details  of  this 
numerical  simulation  below. 

5.2.1.  The  environmental  condition 

The  atmospheric  condition  assumed  for  our  numerical  simulation  of 
melting  layer  is  that  the  temperature  field  satisfies  the  dry  adibiatic  lapse 
rate  from  the  surface  to  the  5°C  level  and  from  the  5^0  level  to  the  0°C 
level,  the  temperature  field  satisfies  the  moist  adiabatic  lapse  rate  with 
100%  RH  with  respect  to  water.  We  selected  a  case  for  surface  pressure  p$  = 
1000  mb  and  temperature  T$  =  10°C. 

5.2.2.  The  computation 

The  generalized  expressions  of  the  cnange  of  the  size  (eqs.  (5)  and  (9)) 
and  the  fall  velocity  (eqs.  (6)  and  (10))  of  snowflakes  or  ice  crystal  under 
constant  temperatures  are  used.  Under  the  atmospheric  condition  we  have 
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chosen,  we  assume  that  within  a  thin  layer  with  thickness  Az  «  25  on  the 
temperature  may  be  regarded  as  having  a  constant  value  ET.  Then,  when  a 
snowflake  or  Ice  crystal  reaches  the  top  of  this  thin  layer  with  the  diameter 
dT  and  fall  velocity  vT,  It  melts  and  falls  under  this  temperature.  There¬ 
fore,  based  on  equations  for  d  and  v,  we  may  obtain  the  initial  size  dQ*  and 
time  t*,  such  that  a  snowflake  or  Ice  crystal  with  Initial  diameter  dQ* 
falling  Into  this  chamber  with  constant  ST  temperature  and  traveling  t*  for 
a  period  of  time  will  reach  a  level  where  they  are  now  dT  and  Vy.  Then,  we 
may  integrate  (6)  or  (8)  to  obtain  the  time  interval  At  required  for  the  snow¬ 
flake  or  ice  crystal  to  travel  through  this  thin  layer  of  Az  *  25  cm,  and 
eventually  from  eqs.  (5),  (6),  (9)  and  (10)  to  obtain  the  size  and  fall 
velocity  when  it  reaches  the  bottom  of  this  thin  layer. 

Based  on  the  calculation  described  above,  we  may  then  follow  the  change 
of  snowflakes  or  ice  crystals  sequentially  downward  until  it  melts  completely 
to  become  droplets.  In  our  calculation,  we  assume  that  the  mass  is  conserved, 
therefore 


Final  density  ,  ^ Initial  diameter 
Initial  density  final  diameter 

As  the  final  density  of  water  droplet  is  1  g/cm9,  while  the  Initial  density 
of  snowflakes  or  ice  crystal  Is  taken  as  0.036  g/cm9.  Hence, 

Final  diameter  «  o  33 
Initial  diameter 

Figs.  40  and  41  describe  the  change  of  the  diameter  of  snowflakes  and 
Ice  crystals  as  a  function  of  altitude  and  Figs.  42  and  43  show  the  fall 
velocity  variation  with  the  altitude.  It  should  be  pointed  out  that  our  fall 
velocity  profiles  are  quite  different  from  that  Austin  and  Bemls  (1960)  guessed. 
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Fig.  40.  Computed  snowflake  dlaroter  d  vs  distance  from  the  0°C  Isotherm  iz 


Fig.  42. 


Instead  of  the  rapid  variation  toward  later  stages  of  melting  In  the  work  of 
Austin  and  Bemls,  the  fall  velocities  we  estimated  vary  with  the  lowered  alti¬ 
tude  In  the  entire  range. 

5.3.  Discussion 

Fig.  44  shows  the  factor  G  computed  In  comparison  with  the  work  by 
Langleben  and  Gunn  with  the  determined  value  of  j  being  32.2.  Fig.  45  de¬ 
scribes  the  radar  reflectivity  computed  from  (67)  with  rf  *  5  mm  and  normalized 
to  that  at  the  0°C  Isotherm.  In  the  computation  It  is  assumed  that  snowflakes 
are  spheres  with  10  mm  diameter  and  obey  the  size  and  the  fall  velocity  the 
regression  equations  give.  There  Is  no  aggregation  taking  place  slmulatane- 
ously.  Also  shown  Is  the  result  with  |K|2  value  for  homogeneous  mixture  of 
water  and  Ice  with  all  other  conditions  being  the  sane.  From  the  figure. 

It  Is  clear  that  the  homogeneous  mixture  treatment  gives  a  much  lower  value 
of  radar  reflectivity  maximum  compared  with  that  of  our  model.  Secondly,  the 
position  of  the  maxlimsn  reflectivity  In  our  model  Is  only  about  40  m  below 
the  0°C  Isotherm  with  the  maximum  4  times  larger  than  that  at  0°C  level.  The 
maximum  reflectivity  position  Austin  and  Bemls  obtained  in  their  model  at 
about  300  m  below  the  0°C  Isotherm  Is  due  to  the  combined  effect  of  several 
factors,  such  as  (1)  rather  sudden  Increase  in  the  assumed  fall  velocity  to¬ 
ward  the  end  of  melting,  (11)  the  gradual  ! K| 2  Increase  due  to  the  homogeneous 
mixture  computation  and  (111)  the  continuing  effect  of  aggregation. 

As  mentioned  In  the  Introduction,  the  radar  reflectivity  in  the  center  of 
the  bright  band  normalized  with  respect  to  that  at  0°C  Isotherm  Is  normally 
between  15  and  30  with  extremes  as  high  as  100  and  as  low  as  2  or  3.  Since 
the  homogeneous  mixture  model  is  Incompatible  with  the  computed  results  of 
Langleben  and  Gunn,  we  compare  our  spherical  snowflake  model  without  aggre¬ 
gation  with  these  observations. 
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Fig.  44.  Variation  of  the  factor  G  with  the  snowflake  diameter  d 


Fig.  45.  The  computed  radar  reflectivity  of  snowflakes  normalized  to  that  at  the  0°C 

Isotherm  n/n*  plotted  as  a  function  of  distance  (temperature)  from  the  Isotherm 
Az  (AT).  Snowflakes  are  assumed  to  be  spherical  with  a  10  mm  diameter. 


There  are  at  least  three  more  factors  to  be  considered  starting  from 
our  model  in  order  to  account  for  the  natural  variation  of  the  bright  band 
conditions,  (i)  According  to  Atlas  et  al.  (1953),  a  partially  melted  ice 
spheroid  gives  a  considerably  stronger  backscattering  compared  with  a  par¬ 
tially  melted  ice  sphere  of  the  same  volume.  It  is  possible  to  expect  a 
similar  effect  with  snowflakes.  Aggregation  of  snowflakes  is  expected  to 
take  place  more  easily  among  partially  melted  snowflakes  leading  to  (11) 
shape  development  and  (ill)  size  increase.  (11)  and  (iii)  are  expected  to 
continue  for  some  time  after  the  spherical  snowflake  reflectivity  has  reached 
its  maximum.  If  three  spherical  snowflakes  form  a  linear  aggregate  at  say 
200  m  below  the  0°C  Isotherm,  it  will  increase  the  reflectivity  by  a  factor 
of  3  from  the  developed  shape,  by  a  factor  of  9  from  the  V2  term  and  by  a 
factor  of  1/3  from  the  number  concentration,  provided  that  the  fall  velocity 
remains  the  same.  This  will  lead  to  an  overall  gain  of  reflectivity  by  a 
factor  of  9.  Multiplying  this  factor  to  the  normalized  reflectivity  of  2 
already  existing  at  the  level,  we  can  expect  the  normalized  reflectivity 
reaching  to  18.  It  will  be  Interesting  to  see  if  the  reflectivity  maximum 
is  larger  and  the  position  lowered  when  the  snowfall  intensity  or  precipi¬ 
tation  rate  is  high  and  vice  versa. 


6.  Concluding  Remarks 

The  present  study  Includes  the  development  of  snow  melting  apparatus, 
the  field  measurement,  the  data  analysis,  the  theoretical  development,  and 
interpretation  of  the  bright  band  mechanism  based  on  the  theoretical  and 
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experimental  knowledge  obtained.  Main  features  and  findings  are  as  fol¬ 
lows. 

The  apparatus  development  phase  Is  highlighted  with  the  following  main 
features : 

(I)  An  air  conditioning  system  to  control  the  temperature  and  relative 
humidity  of  air  passing  through  the  wind  tunnel. 

(II)  A  valve  mechanism  to  control  the  air  velocity  In  the  tunnel  while 
maintaining  the  same  temperature  and  humidity  In  It. 

(11)  A  new  vertically  diverging  wind  tunnel  with  proper  design  and  combi¬ 
nation  of  screens  and  a  honeycomb  to  suspend  snowflakes  and  snow  crystals 
while  melting. 

(1v)  A  trap  door  mechanism  for  sample  snow  Introduction  Ir.to  the  tunnel 
without  causing  turbulence. 

(v)  A  method  to  remove  static  electricity  on  the  tunnel  wall. 

The  main  findings  from  snow  melting  measurements  are: 

(I)  The  observed  melting  time  of  snowflakes  and  snow  crystals  Increased 
with  size,  but  rather  gently,  suggesting  Internal  melting. 

(II)  Snow  crystal  disruption  was  observed  and  recognized  as  a  dominant 
factor  during  the  melting. 

(III)  The  vertical  stability  of  suspended  snowflake  and  snow  crystals 
appeared  to  be  dependent  upon  their  size  and  shape.  Small  crystals,  less 
than  1  mm  In  diameter,  were  the  most  stable. 

(1v)  Helicoptering,  or  the  rotation  of  a  crystal  around  Its  vertical  axis, 
was  observed  In  crystals  1-3  mm  In  diameter.  They  developed  orbital  motion 
with  Increasing  radius  during  the  last  portion  of  their  melting  process. 

(v)  Crystals  less  than  1  mm  In  diameter  juaped  violently  at  the  end  of  their 
melting  process.  The  jump  was  upwards  or  to  the  side,  but  seldom  downwards. 
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(vl)  Large  symmetrical  dendrites  appeared  to  melt  from  the  outside  In  forming 
a  drop  in  the  final  stages  of  melting.  However,  they  too  became  unstable  in 
the  final  stages  of  melting.  There  was  no  disruption  of  these  type  crystals. 

Data  analysis  has  produced  sets  of  regression  equations  for  fall  velo¬ 
city  and  size  as  functions  of  time  and  temperature. 

The  theoretical  development  carried  out  features  with  melting  of  non- 
spherical  single  ice  crystals  and  spherical  snowflakes. 

A  simple  model  of  melting  spherical  snowflake  has  been  developed  in  which 
the  snowflake  is  represented  by  melting  ice  particles  fixed  with  the  space 
coordinate  in  the  snowflake  and  the  snowflake  Is  assumed  to  melt  from  outside 
to  shrink  as  well  as  from  inside  to  incorporate  the  effect  of  water  skin  for 
the  microwave  backscattering.  Comparison  of  this  model  under  the  changes  of 
fall  velocity  and  diameter  predicted  by  the  regression  equations  has  revealed 
that  shape  and  aggregation,  i.e.  size,  of  snowflakes  during  the  melting  stage 
appear  to  account  for  the  observed  increase  of  the  radar  reflectivity  at  the 
center  of  the  bright  band  relative  to  that  at  the  0°C  isotherm  and  the  alti¬ 
tude.  From  this  consideration,  a  relationship  is  predicted  between  the  height 
of  the  center  of  the  bright  band  and  the  snowfall  intensity,  i.e.,  the  lower 
the  center  altitude,  the  stronger  the  snowfall  intensity. 
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